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This paper introduces and studies a field theoretic analogue of the Clebsch 
jy^ ■ variational principle of classical mechanics. This principle yields an alternative 

I ■ derivation of the covariant Euler-Poincare equations that naturally includes co- 

variant Clebsch variables via multisymplectic momentum maps. In the case of 
f— i ' diffeomorphism groups, this approach gives a new interpretation of recently de- 

Q_r rived singular peakon solutions of Diff (M)-strand equations, and allows for the 

^ fH construction of singular solutions (such as filaments or sheets) for a more general 

class of equations, called covariant EPDiff equations. The relation between the 
covariant Clebsch principle and other variational principles arising in mechanics 
i— — i. and field theories, such as Hamilton-Pontryagin principles, is explained through 

the introduction of a class of covariant Pontryagin variational principles in field 
J> ■ theories. 
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1 Introduction 

Variational principles are of central importance in mechanics and field theories, as il- 
lustrated by Hamilton's principle of critical action, stating that the Euler-Lagrange 
equations are obtained by computing the critical curves, or the critical sections, of the 
action functional. On the Hamiltonian side, the canonical Hamilton equations of clas- 
sical mechanics can also be obtained via a variational principle, the so called Hamilton 
phase space principle. This principle extends to the field theoretic context in a natural 
way since the multisymplectic Hamilton equations can be characterized as the critical 
sections of an action functional defined with the help of the Hamiltonian density and 
the Cartan form. 

Variational principles also arise naturally in optimal control theory, via the Pontrya- 
gin maximum principle. Under some assumptions, one can characterize the extremals 
of optimal control problems as the critical points of a functional obtained from the La- 
grange multipliers theorem, e.g. [2]. Such a functional is obtained by adding to the cost 
function of the problem, the control equation constraint multiplied by a costate vari- 
able. This variational principle generalizes naturally to the case of nonlinear optimal 
control problems whose control function takes value in a fiber bundle, e.g. [29]. It will 
be referred to as the Pontryagin variational principle. 

Based on earlier works on geometric optimal control theory such as, e.g., [9], [4], 
[3], [6], [8], [30], [17], a particular class of optimal control problems and variational 
principles associated to group actions was introduced and studied in [24], under the 
name of Clebsch optimal control problem. The main feature of this class of problems 
is that the control function takes values in the Lie algebra of the acting Lie group 
and that the constraint ODE consists in imposing the state of the system to evolve 
according to the infinitesimal action of the control function. The variational principle 
associated to this class of problems, called Clebsch variational principle, provides an 
alternative derivation of the Euler-Poincare equations of classical mechanics based on 
Clebsch variables (as defined by [37]). Moreover, this approach unifies and generalizes a 
wide range of examples appearing in the literature such as the symmetric formulation of 
iV-dimensional rigid body, optimal control for ideal flow using the back-to-labels map, 
double bracket equations associated to symmetric spaces, and optimal control problems 
in imaging science. 
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In a similar way with the Pontryagin variational principle arising in optimal control 
problems, the Hamilton- Pontryagin principle of classical mechanics is a modification 
of Hamilton's principle that includes the Legendre transform as a constraint in the 
action functional, see [39]. This principle combines variables from both the Lagrangian 
and Hamiltonian formulations and is especially appropriate to deal with degenerate 
Lagrangians. Using similar ideas, a field theoretic version of this principle has been 
recently formulated in [40]. 

The first goal of the present paper is to introduce and study a new class of varia- 
tional principles for field theories, called covariant Clebsch variational principle and to 
illustrate the applicability of this principle to a wide range of examples both in the finite 
or infinite dimensional situation. A main example, inspired by the recent work [32], 
being given by the construction of singular solutions for the covariant Euler-Poincare 
equations on diffeomorphism groups. The second goal of the paper is to present the 
unifying geometric framework underlying all the variational principles mentioned above 
(Hamilton phase space principle, Hamilton-Pontryagin principle, Pontryagin principle, 
Clebsch principle) both in the context of classical mechanics and for field theories. 

Plan of the paper. In Section 2 we present some background material needed for the 
rest of the paper. We first review the definition and the main properties of the Clebsch 
variational principle in mechanics. In particular, we explain how the Euler-Poincare 
equations are obtained from the stationarity conditions and mention the occurrence of 
Clebsch variables and collective Hamiltonian. Then we present basic facts about classical 
field theory, such as jet bundles, covariant Euler-Lagrange and Hamilton equations, and 
covariant momentum maps. Finally, we recall the process of covariant Euler-Poincare 
reduction on a trivial principal bundle and relate the spacetime decomposition of these 
equations with the G-strand equations of [32]. In Section 3 we introduce the covariant 
Clebsch variational principle and study its main properties. In particular, we present 
the stationarity conditions, explain their covariant Hamiltonian formulation through 
the introduction of a collective Hamiltonian via the covariant momentum map, and 
show that they yield the covariant Euler-Poincare equations. In Section 4 we apply 
the covariant Clebsch variational principle to various particular cases of special interests 
such as the case of a trivial fiber bundle, the case of group translations, and the case 
of a linear action. In the case of diffeomorphism groups, this approach gives a new 
interpretation of recently derived singular peakon solutions of Diff (M)-strand equations, 
and allows for the construction of singular solutions for a more general class of equations, 
called covariant EPDiff equations. We also consider the covariant version of the Euler 
equations for perfect fluids, called covariant EPDiff,,^ equations and provide covariant 
Clebsch variables for these equations. Finally, in Section 5 we present a class of covariant 
variational principles for field theories that underlies the covariant Clebsch variational 
principle and other known principles such as the Hamilton phase space principle and 
the Hamilton-Pontryagin principle, both in their classical and covariant form. 
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2 Preliminaries 



In this section we review some needed facts concerning the Clebsch variational principle 
in mechanics and the multisymplectic geometry of classical fields. 



2.1 The Clebsch variational principle in mechanics 

Recently, a particular class of geometric optimal control problems associated to Lie 
group actions has been considered and studied in [24], see also [17], in relation with 
variational principles and Clebsch variables in mechanics. The main feature of this class 
of problems is that the control function takes values in the Lie algebra of the acting 
Lie group and that the constraint ODE consists in imposing the state of the system to 
evolve according to the infinitesimal action of the control function. It was shown that the 
necessary condition for an optimal solution, obtained via the Pontryagin maximum prin- 
ciple, naturally encodes a Clebsch representation of the momentum variable conjugate 
to the control. The associated variational principle is therefore naturally referred to as 
the Clebsch variational principle. It is reminiscent of the Hamilton-Pontryagin principle 
which incorporates, instead of a Clebsch representation, the Legendre transformation, 
and enables one to treat degenerate Lagrangian systems, see [39] and references therein. 
Both variational principles can be seen as particular instances of the variational principle 
for optimal control problems provided in [29] in connection with Pontryagin maximum 
principle. This latter variational approach will be adapted to the covariant field theoretic 
context in the last section of the present paper. 

Various optimal control formulations of mechanical systems were shown in [24] to 
arise as particular instances of Clebsch optimal control problems, obtained by an ap- 
propriate choice of the Lie group action on the state manifold. Such examples include 
various double bracket equations and systems on Lie algebras and symmetric spaces ([9], 
[4], [3]); optimal control problems for fluids ([5], [30], [17]); the symmetric representa- 
tion of the rigid body and its generalizations ([4], [6], [8], [7]). A particularly relevant 
Clebsch representation for the present paper is the one provided by the singular solution 
of the Camassa-Holm equations and its higher dimensional extensions ([34], [17]). We 
now recall the definition of the Clebsch optimal control problem. 



The Clebsch optimal control problem. Let $ : G x Q — )■ Q, (g, q) i— > $ g (q) be 
a left action of a Lie group G on a manifold Q. Let us denote by q the Lie algebra of 
the Lie group G and by X(Q) the space of all smooth vector fields on the manifold Q. 
Given a Lie algebra element £ G g, the associated infinitesimal generator £q G X(Q) is 
defined by 

$cx P (to(<?) ? q^Q, 

t=o 

where exp : q — > G is the exponential map of G. Given a cost function £ : q — > R, a real 
number T > 0, and two elements qo, qr G Q belonging to the same G-orbit, the Clebsch 
optimal control problem associated to the action $ consists in minimizing the integral 

1 mt))dt (2.1) 



fc (ff): dt 
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subject to the following conditions: 



(A) q(t)=mQ(qm 

(B) q(0) = q and q(T) = q T . 

From the Pontryagin maximum principle, it follows that an extremal solution (q(t), £(£)) 
of (2.1) is necessarily a solution of the Clebsch variational principle 

S f 1 \(a,q-t Q (q))+£(O)dt = (2.2) 
Jo 

over curves £(£) G and a(t) G T*Q, where a(t) is a curve covering g(t) G Q with fixed 
endpoints go an d qr- This principle yields the stationarity conditions 

(see [24]), where J : T*Q — > q* is the momentum map associated to the cotangent lifted 
action of G on T*Q, and £t*q G X(T*Q) is the infinitesimal generator of this cotangent 
lifted action. Recall that J is given by (J(a q ),£) := (ct 9 ,£Q(g)), for all £ G g, a q G T*Q. 
Recall also that the cotangent lifted action of $, denoted by $ T *, is defined by 

<S> T * : G xT*Q ^T*Q, $f (a,) := T*$ 9 -i(a,), for all a q eT*Q, 

where T*$ 5 -i : T*Q -»■ T*g is the dual of the tangent map T$ g -i : TQ TQ. 



Recovering the Euler-Poincare equations. A solution of the second equation in 
(2.3) is necessarily of the form a(t) = $J^(a(0)), where ^(t)p(t) -1 = Using the first 
equation in (2.3), we obtain that the control satisfies the Euler-Poincare equations 

d 5£ „ 51 

where 5£/5£ G g* denotes the functional derivative of £ and (ad^ /x, r/) = (/x, [£,77]), with 
£,77 G and /x G g*, is the infinitesimal coadjoint action. Recall that the Euler-Poincare 
equations can be obtained by Lagrangian reduction of the Euler-Lagrange equations 
associated to a Lagrangian L : TG — > M defined on the tangent bundle of G and G- 
invariant under the tangent lift of right translations. A curve g(t) G G is a solution of the 
Euler-Lagrange equations for L = L(g, g) if and only if the curve £(£) = g{t)g{t)~ l G is 
a solution of the Euler-Poincare equations for the reduced Lagrangian I := L\ s : — > M. 
The same procedure on the Hamiltonian side leads to the so called Lie-Poisson equations 

ji + adk 7x = (2.4) 

S/j, 

associated to a given Hamiltonian h : g* — > R. These equations can be obtained by 
Poisson reduction of canonical Hamilton's equations for the right invariant Hamiltonian 
H on T*G with h = H\ s . If £ \— > J| is a diffeomorphism, then Euler-Poincare and 
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Lie-Poisson equations are equivalent, since one can pass from one to the other by the 
Legendre transformation 

5£ 

h(M) = M-t(t), H = jz. (2.5) 
Let us recall that the equations (2.4) are Hamiltonian relative to the Lie-Poisson bracket 

fi e q* (2.6) 

on functions /, g : g* — > R. We refer to [36] for a detailed account about Euler-Poincare 
and Lie-Poisson reductions. 

Collective Hamiltonian and Clebsch variables. The Pontryagin Hamiltonian as- 
sociated to the Clebsch optimal control problem (2.1) reads 

H(a q ,£) = (a, t Q (q)) -£(£). (2.7) 

Assume that £ i— >■ J| is a diffeomorphism and denote by £*(cK g ) the optimal control 
uniquely determined by the condition S£/S^ = 3(a q ). The optimal Hamiltonian of the 
problem, defined by H(a q ) := H(a q ,£*(a)), is thus given by 

H(a q ) = (|,e) -^(0 = h (£) = h(J(a q )), (2.8) 

where h : g* — > R is the Hamiltonian associated to £ via the Legendre transformation 
£ i — y S£/S^, see (2.5). Therefore, in the case of Clebsch optimal control problems, the 
optimal Hamiltonian coincides with the so called collective Hamiltonian associated to 
Lie-Poisson Hamiltonian h and the momentum map J : T*Q — > q*, (see [28], [36]). 

We now recall from [37] the abstract notion of Clebsch variables. Let P be a Poisson 
manifold, thought of as the phase space of a physical system. By Clebsch variables (or 
symplectic variables) for this system, we mean a symplectic manifold R and a Poisson 
map 

* : R -> P. (2.9) 

Classical examples of such Clebsch variables are the Cayley-Klein parameters for the free 
rigid body, the Kustaanheimo-Stiefel coordinates in quantum mechanics, and, of course, 
the classical Clebsch variables in various fluid dynamical models (homogeneous incom- 
pressible, isentropic, MHD, etc), see [33]. Since the momentum map J : T*Q — > q* is a 
Poisson map with respect to the canonical symplectic structure on T*Q and the (non- 
canonical) Lie-Poisson structure (2.6) on q*, the optimal control condition 5£/5£ = J(ce q ), 
naturally defines a Clebsch representation for the momentum 5£/5£. More precisely, 
the noncanonical Hamiltonian motion of the control is found to arise from canonical 
Hamilton's equations for the collective Hamiltonian, via a momentum map. Such an 
observation is particularly relevant for the dynamics of the singular solutions of the 



S_l Sg_ 
<5/x' 5jjL 
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Camassa-Holm equations and their generalizations (called EPDiff), which arise from 
the momentum map 

(Q, P) G T* Emb(S, M) h+ m = / P(s)5(x - Q{s))ds G £(M)*, 

Js 

([34]) as we will recall later. This expression of m in terms of (Q, P) is a Clebsch repre- 
sentation which yields a singular solution of the EPDiff equations in terms of canonical 
variables that obey canonical Hamilton's equations on the cotangent bundle of the space 
Emb(S, M) of all embeddings Q : S ->■ M. 

2.2 The geometry of jet bundles and classical field theory 

In this section, we recall some basic facts about the geometry of classical field theory. 
This includes the notion of jet bundles, dual jet bundles, canonical multisymplectic 
forms, covariant momentum maps, covariant Euler-Lagrange and Hamilton equations. 
For detailed informations, the reader can consult e.g. [26], [27], [18], [11] and references 
therein. 

Jet bundles. Let ttx.y '■ Y — > X be a fiber bundle which we call the covariant con- 
figuration bundle. The fiber at x G X is denoted Y x = n x l Y (x). The physical fields are 
sections of this bundle, that is, smooth maps <fi : X — > Y such that ir o = idx- We 
suppose that dim X = n + 1 and dim Y = n + N + 1. Coordinates on X are denoted x^, 
\i = 0, ...,n, and fiber coordinates on Y are denoted y A , A — 1, ...,N, so that a section 
reads locally ^(x^) = (x M , y A (x^)). 

The first jet first bundle tty.j^y '■ J X Y — > Y is the affine bundle over Y whose fiber 
at y G Y x is 

J]Y := {j y G L(T X X, TyY) \ Ttt x ,y ° l y = id Tx x} , 

where L(T x X,T y Y) denotes the space of linear maps T X X — > T y Y . It is the covariant 
analogue of the tangent bundle in classical mechanics. Given a section of kx,y '■ Y — > 
X, the map x i— > j x 0(x) := T x (p defines a section of the fiber bundle ttx,j 1 y '■ J l Y — > X. 
This section is called the first jet prolongation of <fi and reads locally 

x^ 1 ^ (x^,y A (x^ J d„y A (x^)). 

Dual jet bundles. The field-theoretic analogue of the cotangent bundle or phase 
space of classical mechanics, is given by the dual jet bundle tiy,j^y* '■ J X Y* — > Y, defined 
as the vector bundle over Y whose fiber at y G Y x is the space 

J^ = Aff(J y 1 F,Ar 1 X) (2.10) 

of all affine maps JyY — > A X +1 X, where A n+1 X is vector bundle of (n + 1) forms on 
X, supposed to be an orientable manifold. Any affine map from JyY to A X +1 X can be 
locally written as 

where d n+1 x = dx°Adx 1 A...Adx n , so that coordinates on J X Y* are given by (x M , y A ,p^, vr). 
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The dual jet bundle ti y ,j^y* '■ J X Y* — > Y is canonically isomorphic to the vector 
bundle tty,z '■ Z — > Y, whose fiber at y G Y is given by 

Z y = {z y e Ky +1 Y | \ v \ w z y = 0, f or all v , w G V y Y } , (2.11) 

where denotes the interior multiplication by v and V y Y := ker(T y 7Tx,Y) C T y y is the 
vertical subspace. The canonical isomorphism is provided by the vector bundle map 

$ : Z ->■ J 1 ! 7 '*, {§{z y ),7 y ) := <y*z, for all z y G Z y and 7 y G J^F. 

Locally the inverse of this isomorphism reads 

+ P^fy (8) (fy A ) d n+1 x ^ vrrf n+1 + p A dy A A d%, 

where := \Q^d n+1 x. 

The restricted dual jet bundle hyti '■ n — )■ F is, by definition, the vector bundle over 
Y whose fiber at y is 

n, := L(L(T x X,V y Y), A n x +1 X). (2.12) 

The total space II can be seen as the quotient of JyY* by constant affine transformations 
along the fibers of vry jiy, so that coordinates on II are given by (x M , y A ,p'X) and the 
quotient map ji : J X F* — > II reads locally 

Note that fi : J l Y* -» n is a line bundle, the fiber being isomorphic to A™ +1 X. 

Canonical multisymplectic Forms. Analogous to the canonical forms Q can = PAdq A 
and fl C an — dq A A dp a on a cotangent bundle, there are canonical forms on the dual jet 
bundle J X Y* ~ Z . The canonical (n + l)-form G ft n+1 (Z) is defined by 

e(z)(u ,...,u n ) := z(T-k y ,jiy*{u q ), ...,T-K Y ^ Y *(u n )) = (7i Y JlYi z)(u ,...,u n ) 

and the canonical multisymplectic form G fl n+2 (Z) is defined by 

ft := -d0. 

These forms are locally given by 

= PA^y A A d n x fl + ird n+1 x and ft = A/ A A dp^ A d% - dix A rf n+1 x. (2.13) 

Lagrangian formalism. Assume that X is an oriented manifold. A Lagrangian den- 
sity is a smooth map £ : J X F — )■ A n+1 X covering the identity on X. In local coordinates, 
we can write 

C(x^y A ,v A ) = L(x^y A ,v A )d n+1 x. 

Let £/ C X be an open subset whose closure U is compact. Throughout the paper, 
we will assume that the boundary dU is regular enough to enable the use of Stokes' 
theorem. Recall that a local section : U C X — > Y of itx,y is, by definition, smooth 
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if for every point x G U there is an open neighborhood U x of x and a smooth section 
<Px '■ U x — > Y extending 0. A critical section of the action functional defined by £ is a 
smooth local section 6 : U Cl->y that satisfies 



d_ 

de 



£(jV £ ) = 0, (2.14) 



--0-JU 



for all smooth variations <p £ : U C X — > Y such that </>o = and 4> £ \du — 4>\du- 

It is well-known that a section is critical if and only if it verifies the covariant Euler- 
Lagrange equations 

in local coordinates. 

Hamiltonian formalism. There are two equivalent field-theoretic analogues of the 
Hamiltonian function of classical mechanics: the Hamiltonian density and the Hamilto- 
nian section. 

For the approach using the Hamiltonian density, we follow the formulation given in 
[11] (see in particular the Appendix), which is independent of the choice of an Ehresmann 
connection on tix,y '■ Y — > X or a volume form on X. By definition, a Hamiltonian 
density is a smooth map H : J X Y* — > A n+l X covering the identity on X such that 

i x (dH + Q) = 0, for all X G X V (J 1 Y*), (2.15) 

where X v \J l Y*) = {X G 3L(J l Y*) \ Tfi o X = 0} is the space of all //-vertical vector 
fields on the dual jet bundle. From the condition (2.15) we obtain that a Hamiltonian 
density is locally given by 

U{x^y A i P %Tt) = (n + h loc (x^y A ,p» A ))d n+1 x. 

A Hamiltonian section is, by definition (see [26]), a smooth section h : n — > ^Y* of 
the line bundle \i : J l Y* — > n. It is therefore locally written as 

h{x^y A y A ) = (x^y A ,p^-h loc (x^y A ,p'X)). 

As the local expressions suggest, there is a bijective correspondence between Hamil- 
tonian densities and Hamiltonian sections. This correspondence being given by 

imh = n^(0). 

The Cartan (n+ l)-form and (n + 2) forms associated to a Hamiltonian section h are 
defined by 

6 h := h*0 G fi n+1 (n) and fi h := h*tl G n n+2 (U) 
and we have the relations 

fj,*e h = <d-H and = Vt + dU. 
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On the Hamiltonian side, the physical fields of the theory are represented by local 
smooth sections u : U C X — > II of the restricted dual jet bundle 7ix,u '■ n — > X. The 
covariant analogue of the classical Hamilton's phase space variational principle reads 



d_ 

de 



e 



u* e e h = 0, (2.16) 



for all smooth variations u e : U C X — > U of a given section u : U — > II, such that 
Wo = w an d </> e |d[/ = </>|dtr, where £ := 7i Y ,n ° uj £ : U Cl->7. A section u; is critical if 
and only if it verifies the covariant Hamilton equations 

(j/(iyf2h) = 0, for every 7Tx,n-vertical vector fields V on IT. 

In local coordinates, these equations take the standard form 



du A _ dh loc dtJX _ dh 

dx^ du A ' dx^ du 



TTV = -VfM- (2-17) 



The Legendre transformations. Given a Lagrangian density C : J X Y — > A n+1 X, 
the associated covariant Legendre transformation is defined by 



F£ : J l Y -> J 1 ^, (F£( 7 ), 7 ') := £(t) + ^ 



£( T + £(7-7')), 



£=0 



where 7, 7' G J* F for some y EY . The reduced covariant Legendre transformation is, 
by definition, 

FCi^F^n, ¥C:=fio¥C. 
In local coordinates, they read 

^ O = ^ j|> * - j|<) and O = ^ j|) • 

The Lagrangian density is said to be hyperregular when F£ : J 1 !" — > U is a diffeo- 
morphism. In this case, it is possible to define the Hamiltonian section associated to L 
as 

h c := F£oF£ _1 : n -»■ ^F*. 

It is locally given by 

hc{x^y A M = (x»,y A ,p%L(x»,y A ,v A )-p A v A ), 

where, by the hyperregularity hypothesis, v A is uniquely determined from the other 
variables by the condition Jj^- = p A . The associated Hamiltonian density %c : J X Y* — > 
A n+1 X is given locally by the expression 

n c (x^y A ,p%7r) = (tt + p A v A - L(x»,y A ,v A ))d n+1 x. 

In the hyperregular case, one checks that the covariant Euler-Lagrange and covariant 
Hamilton equations are equivalent. 
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Covariant momentum maps. Let Q be a Lie group acting by bundle automorphisms 
on the left on the configuration bundle ttx,y '■ Y — > X. Let us denote by ip : Q x Y — > Y 
this action and by ip x : Q x X — > X the induced action on X. This action naturally 
induces ^-actions on J l Y and ,PY* . They are respectively given by 

tfY : jiy jiy; ^( 7y ) := Tip g o 7, o ^ 

and 

where we identified J X Y* with Z. The lifted action on the dual jet bundle is a special 
covariant canonical transformation in the sense that its preserves the canonical (n + In- 
form and, therefore, the multisymplectic form Q: 

(i/>f Y *)*e = 6 and (i/)f Y ySl = n, for all g G Q. 

Recall ([27]) that the covariant momentum map associated to the lifted action ty J g Y * 
is given by 

J : J x y* L(g, A n J 1 y*), J(z)(0 = n* YtZ is Y z, (2.18) 

where £y G X(Y) is the infinitesimal generator associated to the action of Q. It is locally 
given by 

A*)(0 = (Pa^ + ^ M M% - ^ V A <f*~ V> (2.19) 
where d n ~ 1 x fll/ := ig„ig d n+1 x and ) is the local expression of £y. The map J 

verifies the momentum map condition 

where £jiy* is the infinitesimal generator of the lifted action on ^Y* and J^ G f2 n+1 ( J x y*) 
is the (n + l)-form defined by J^(^) := J(z)(£). 

2.3 Covariant Euler-Poincare reduction and G-strands 

In this paragraph, we first recall the process of covariant Euler-Poincare reduction on a 
trivial principal bundle. It extends the Euler-Poincare theory recalled above in §2.1, to 
the field theoretic context. Covariant Euler-Poincare reduction for principal bundle field 
theories has been developed by [14], and was then extended to more general situations in 
[13] and [19]. In the second part of this paragraph, we relate the spacetime decomposition 
of these equations with the continuum spin chain or G-strand equations considered in 
[31], [32]. 

Covariant Euler-Poincare reduction. Suppose that the covariant configuration 
bundle is a trivial right principal bundle ttx,p '■ P = X x G — > X. Since the bun- 
dle is trivial, the fiber J 1 P( 2;)9 ) of the first jet bundle can be canonically identified with 
the vector space L(T X X, T g G). Sections of this bundle are identified with smooth maps 
g : X — > G, via the relation <j)(x) = (x,g(x)). 

Let C : J X P — > A n+1 X be a Lagrangian density supposed to be right G- invariant, that 
is C(TR h o V( x ,g)) = £(u( Xjff )), for all h G G, v^ x>g ) G L(T x X,T g G), where TR h : TG -> 
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TG denotes the tangent map to right translation by h on G. Since the quotient 
vector bundle J x PjG — > X can be identified with the vector bundle L(TX, g) — > X, 
the map C induces a reduced Lagrangian density £ : L(TX, g) — > A n+1 X defined by 
- 1 ° ^(g,^)) — £-{v(x,g))- The relation between the field g : X — y G of the theory 
and its reduced expression cr is a = TR g -i o dg, where dg : TX — >■ TG is the tangent 
map to g. We will use the simpler notation a = dgg -1 . Hamilton's variational principle 
(2.14) for C induces the following constrained variational principle for £: 



5 / £(a) = 0, for all 5a = d( + ad f a, (2.20) 
Ju 

where a : U C X — > L(TX,q) is a local section and ( : U C X — )■ g is an arbitrary 
smooth map such that (\ d jy = 0. The stationarity condition yields the covariant Euler- 
Poincare equations 

5£ 5£ 

div^ + ad;-^ = 0, (2.21) 
5a 5a 

where y- € L(L(TX, g), A n+1 X) is the functional derivative of £. In addition to the 
covariant Euler-Poincare equation, the reduced field a also satisfies the zero curvature 
condition 

da=[a,a]. (2.22) 

Conversely, in order to reconstruct from (2.21) a solution g : U C X — > G of the 
covariant Euler-Lagrange equations for C, one has to assume the zero curvature condition 
(2.22). 

We refer to [11] for the detailed derivation of these equations together with the full 
treatment on arbitrary principal bundles. 

Spacetime decomposition and G-strand equations. In the special case when the 
base manifold decomposes as X = K. x M 3 {t,m), we can write the reduced field as 
a(t,m) = u(t,m)dt + 7(t,m), where v(t,m) G g and 7(t, m) G L(T m M, g). Note that 
the relation a = dgg^ 1 between the initial field g and its reduced expression yields the 
equalities 

v = dtgg^ 1 and 7 = d M gg^, 

where dM<? : TM — > TG denotes the partial derivative with respect to the variable 
m G M. By writing explicitly the covariant Euler-Poincare variational principle (2.20) 
in this case, we get 

5 / £(u, 7) = 0, for all 5v = d t ( + v and £7 = dj^C + a d§ 7, (2.23) 
Ju 

where we choose U —]ti, t 2 [xV with V an open subset in M and ( : U — > g is an arbitrary 
map with = 0. The covariant Euler-Poincare (2.21) can thus be equivalently written 

as 

d 5t 5t 51 51 

at 5v 5u 5^ ' 5^ 

and the zero curvature condition becomes 

d t j - d M v = [v, 7] and d A /7=[7,7], (2.25) 
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see [23]. 

In the case M — R, the equations (2.24), together with the variational principle 

(2.23) , coincide with the ones considered in [31] (see §10) in connection with the dynamics 
of a continuum spin chain. These equations, also called G-strands, were further studied in 
[32] for several remarkable choices for the group G, in connection with integrable classical 
chiral models. The case G = Diff(R), where Diff(R) denotes the diffeomorphism group 
of the real line, is especially important. Indeed, [32] show that when the Sobolev H 1 
metric is used in the Lagrangian, the Diff(M) -strands equations admit peakon solutions 
in a similar way as the Camassa-Holm equations. We will recall this fact in detail 
later (in §4.4) and extend the peakon solution to sheets and filaments. Moreover, we 
will provide a covariant momentum map interpretation of these singular solutions, and 
prove that they undergo canonical covariant Hamiltonian dynamics. 

Note that in the case of G-strands, denoting R x M = R x R 3 (t,s), equations 

(2.24) and the first equation in (2.25) reduce to 

d S£ 5C d 81 81 

whereas the second equation in (2.25) vanishes. 

Remark 2.1 (A classical Lagrangian reduction approach) Remarkably, the covari- 
ant Euler-Poincare equations (2.23) can also be obtained via a classical (as opposed to 
covariant) Lagrangian reduction approach in the sense that the field 

g : X = R x M -> G, (t, m) h> g(t, m) 

is interpreted as a curve g : t G R — > g(t) G J^(M , G) in the infinite dimensional 
Lie group J^(M, G) rather than a section (t,m) G R x M — > (t,m, g(t,m)) G R x 
M x G of the principal bundle. This approach (developed in [23]) makes use of the 
process of affine Euler-Poincare reduction for the Lie group J-"(M, G) acting by affine 
representation on the space f2 1 (M, g) 3 7. In this case, the first equation in (2.25) 
becomes an advection equation, while the second equation is a consequence of the initial 
value (zero, for instance) chosen for the variable 7 and can be generalized to include 
nontrivial dynamic of the curvature B = d-y — [7,7]. We refer to [22] for a detailed 
treatment of affine Euler-Poincare reduction, and to [23] for its application to covariant 
Euler-Poincare reduction. 



3 The covariant Clebsch variational principle 

3.1 Covariant formulation of the standard Clebsch principle 

In this section we formulate the classical Clebsch variational principle reviewed in §2.1, 
by using the geometric tools of covariant field theories recalled in §2.2. This will be 
crucial to introduce the covariant Clebsch variational principle in the next paragraph. 
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Recall from §2.1 that, for a given Lie group action $ : G x Q — > Q and a Lagrangian 
I : q — > R, the Clebsch variational principle is 

5 [ T (£(£) + {a, q-£ Q (q)))dt = (3.1) 
Jo 

over the space of smooth curves (£, a) : [0, T] — > X T*Q, such that g(0) = go and 
g(T) = q?, where q(t) e Q is the curve induced by a(t) G T*Q on Q. We now rewrite 
the integrand in a covariant form suitable for the extension from R to an arbitrary base 
manifold X. Defining the function e(a q ,t,£) := (J(a g ),£) — we can write 

(£(0 + K g - £ Q {q)))dt = m) - (JK), £> + K, <7» 

= — e(a , t,£)dt + (a, g) c/t 
= fe*9(f), 

where the map ip, defined by ip(t) := (az(t),t,£(t)) is interpreted as a section of the 
trivial bundle 

the function e, defined by e(a q ,t, £) := (a 9 , i, — e(i, a g , £), £) is interpreted as a section 
of the trivial bundle 

T*(Q x R) x T*Q x R x , ( a „ t, vr, ^ (a,, t, £); 

and the one-form 9 on T*(Q x M) x g is defined by 0(a o , t, 7r) = ndt + 6 can (a g ), where 
©can £ fi 1 (T*Q) is the canonical one-form on T*Q. 

At this point one recognizes that for the covariant configuration bundle of classical 
mechanics, namely 

Y = R x Q -> R, (t,q)^t, 

we have (from (2.10) (2.12)) 

J 1 ^ = T*(Q x R) 9 (a 9 ,t,7r) and n = T*Q xl9 (« ,t). 

We now incorporate the Lie algebra in each F-fiber by considering the bundles 

J 1 ^* : = J 1 !-* Xy and II B := IT x y (3.2) 

over Y = Q x R, so that ip and e can be seen, respectively, as sections of the bundles 
n g — > R and J l Y* — > R. The variational principle (3.1) can thus be rewritten as 

5 f ^* e *e = 0. 
Jo 

It is this geometric formulation that we will use to define the covariant Clebsch varia- 
tional principle in the next paragraph. 

Note also that the action of G on Q naturally induces an action by bundle automor- 
phisms on Y = Q x R obtained by acting trivially on R. The covariant momentum map 
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J associated to this action recovers the usual momentum map J on T*Q since, from 
formula (2.19), we get 

J(a q ,t,ir) = 3{a q ). 

Therefore, the section e can be expressed in terms of the covariant momentum map as 
follows 

e(a„t,0 := (a q ,t, -e(t, a„ e (^"^0 = (J(<Xq,t,n),£) - £{£)■ (3.3) 

3.2 The covariant Clebsch variational principle and its proper- 
ties 

3.2.1 Definition of the covariant Clebsch variational principle 

In order to define the covariant Clebsch variational principle, we fix a fiber bundle 
tix,y '■ Y — y X and an action ip g : Y — > Y of a Lie group Q. We suppose that tp g covers 
the identity on X, that is hx,y ° — ^x,y- The coordinates of a Lie algebra element 
£ G g will be denoted by £ Q and the coordinates of the associated infinitesimal generator 
£r by (e,£ A ). 

Let L(TX,g) — y X be the vector bundle whose fiber at x consists of linear maps 
a x G L(T X X, g), where g is the Lie algebra of Q, and fix a map £ : L(TX, g) — y A n+1 X 
covering the identity on X. 

Such a map can be obtained by reduction of a (/-invariant Lagrangian density C : 
J X P — y A n+l X, defined on the first jet bundle of the trivial principal bundle P := 
X x Q — y X, as recalled in the last paragraph of §2.2. 

As in (3.2), we need to extend the line bundle \x : J l Y* — y II in order to incorporate 
the vector bundle L(TX,g) —y X. This is done by considering the line bundle 

fi s : J 1 ^ ->■ n s , where J l Y* := ^Y* x Y n XiY L(TX, g) 

II g :=Ux Y ir* X)Y L{TX,g) 

In order to formulate the variational principle, one needs to define a section e of this line 
bundle which appropriately generalizes the expression (3.3). We will use the following 
reformulation of the covariant momentum map of the (^-action, namely, we define the 
map 

J : J'Y* L(g, A n Y), J{z){i) := i^z. (3.4) 

We thus have the relation J(z)({,) = n YiJ i Y *(J(z)(0), where J : J l Y* -> L(g, K n J l Y*) 
is the covariant momentum map as defined in (2.18). 

The construction of the section e will use the following lemma. 

Lemma 3.1 Let Q be a Lie group acting on hx,y '■ Y — X by bundle automorphisms. 
Then, for all a x G L(T X X, g) and z y G Z y ~ JyY*, where nxyiy) = x, we have 
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where, by definition, the wedge product involves a contraction of the Lie algebra indices. 
We thus obtain a well defined bilinear bundle map 

j : J l Y* = J'Y* x Y tt xy L(TX, g) -> J l Y\ ](z y ,a x ) = n* x>Y (a x ) A J(z) 

covering the identity on Y . 

When the Q -action covers the identity on X , this map induces a well defined bilinear 
bundle map 

j : n fl = n x y tt XjY L(TX,q) J l Y* (3.5) 

such that j o ft = j . 

Proof. Since n XtY (a x ) G L(T y Y,g) and J(z) G L(q, A™Y), we have ir XY (a x ) A J(z) G 
A n+1 Y. So, to prove that this is an element in the subspace Z y C A™ +1 y, we need 
to check that it verifies the definition (2.11). This is trivially the the one-form 

7T XY (a x ) annihilates vertical vectors. 

If the action covers the identity on X, the infinitesimal vector field £y is vertical. 
Thus, from the local expression of the momentum map (2.19), it follows that J does 
not depend on the coordinates ir in any local charts. This proves that the expression J 
only depends on the projection \i{z) G II of z. We thus obtain the desired result. ■ 

In order to define the section e, we need to fix an Ehresmann connection on the fiber 
bundle kx,y '■ Y — > X, that is, a smooth subbundle HY of TY called the horizontal 
subbundle such that TY = HY © VY. This is equivalent to the choice of a smooth VY- 
valued one-form 7 G 1 (Y', VY) such that •y\vY = idyy- in coordinates, the evaluation 
of 7 on a tangent vector to Y, namely (v M , v A ) is written as (0, v A + 7^w M ). As shown in 
[26], such a connection defines a linear section s 7 of the line bundle fi : J l Y* — > U given 
by 

s 7 : n = L{L{TX, VY),A n+1 X) -> J X Y\ s^v <g> a ® /x) := 7*a A i v fx, 
where a G VY*, \i G A n+1 X, and v G TX. Locally, the section reads 

s^x^y A M = {x^,y A y A ,p^). 

By assuming that the ^-action covers the identity on X, we can now define the map 
e as 

e(w, a) := (s 7 (w) - a) + 7T* X Y £(a), a) (3.6) 
and we have the following result. 

Lemma 3.2 The map e is a section of the vector bundle 

N : j%* -> n . 

It is locally given by 

e(^, y A ,p\, ol) = (x», y A ,p% -p» A (a A - tf) + l(o$, <), (3.7) 
where I is the local expression of £ defined by £((t") = l(a")d n+1 x. 
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Proof. It suffices to show that fi g o e = idn . We have 



JU B (e(w, a )) = (M s7 M - i( w ? CT ) + 7r x,y^( (T ))> CT ) 
= (u,a), 

where we used fi(s 7 (u)) = u which follows from the fact that s 7 is a section of the line 
bundle \i : J l Y* — > IT, and the equalities /i(j(cu,a)) = 0, /i(7r^ y £(cr))) = which both 
follow from the property /z^* y au+i x = 0. To see that j(u,cr) is the pull back by ttx,y 
of a (n + l)-form on X, we use local coordinates. We have, locally, 

a) = a a A ,7(*)a = c^cfe^ A J{z) a = dx^ A J(z){a ll ) 
= dx» A p v A o*d n x v = p^d n+1 x, 

which yields the desired result. 

The local expression (3.7) follows from the local expression of each term. ■ 

Making use of the above considerations we can now state the main definition of this 
section. 

Definition 3.3 (Covariant Clebsch variational principle) Let tt x ,y : Y X be a 
fiber bundle endowed with an Ehresmann connection 7, let Q be a Lie group acting on Y 
and covering the identity on X, and let I : L(TX, g) — > K n+1 X be a Lagrangian density. 
Define the (n + 1) form 

where e is the section associate to i and 7 defined in (3.6). Let U C X be an open 
subset whose closure U is compact, and let ip : U C X — > U Q be a local smooth section 
°f ^x.rig : n g X. The covariant Clebsch variational principle is 

5 I ^*6 e = 0, (3.8) 
Ju 

for all variations ip E : U C X — > IT x Y L(TX, g) ofip (among smooth sections) such that 
ipQ = ip and (j) s \du — 4>\du> where <p e is the section of ttx,y -Y^X induced from if) e . 

3.2.2 Stationarity conditions 

In this paragraph we compute the equations characterizing the critical points of the 
variational principle (3.8). We present the computations both in global formulation and 
in local coordinates. 

Proposition 3.4 A smooth section ip : U C X — > U B is a critical point of the covariant 
Clebsch variational principle if and only if it satisfies 

ip*(iv^e) = 0, for any iix,u B -vertical vector fields V on U Q , (3.9) 

where Q e := — dQ e . 
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Proof. Since 

ijj e (x) G V^ x )U a := ker(T^, (a; )7rx J n e ) and 5ip\ du = 0, 



e=0 



one may assume without loss of generality that ip e = r] £ o ip, where rj E is the flow of a 
7Tx,n e -vertical vector field V G X^Y*) such that V(^(a;)) = for all x G dU. Let 
us denote by £ya the Lie derivative of a /c-form a. Recall that it verifies the Cartan's 
magic formula £ya = iyda + diy<x It follows that 



0(£ 



4*e e = / ij*£ v ®e 

z=0JU JU 



[ V*(-ivfte + d(i v e e )) 



= - / ^*iyfi e + / ^i^6e = - / ^*iy^ e , 

J(7 ./at/ Ju 

for all V, where we utilized Stokes' theorem and V(ip(x)) = for all x G dU. Thus ^ is 
a critical point if and only if (3.9) holds. ■ 

We now compute the stationary condition in local coordinates. 
Lemma 3.5 The stationary condition (3.9) reads locally 

|g=«(l?-0,.. % = ^f'^W. (3-10) 

where we denoted by ip(x^) = (x^ , ^ A (x^) , ip^x^ 1 ) , a^(x^)) the section ip in local coordi- 
nates. 

Proof. Using the local expression of the multisymplectic form Q, see (2.13), we get 

n e = dy A A df A A d% + (j^dy A + + ^ d <) A ^ 

where e/ oc (x M , y A ,p A , o~") = Pa(o~ a ~ it) ~ K^)- Note that a 7Tx,n B -vertical vector field 



on IT g reads locally V = (V = 0, V A , V%, V«). We thus have 
i v n e = V A df A A d% - V^dy A A d% + (^f ^ + |^V£ 



d&loc t ra \ m+l 



and, therefore, 



Since the vector field V is arbitrary, we get the conditions 



_ dei oc d^ A _ dei oc de 



loc 



g x fi QyA ' Q x fi Qpl^ ' Q a a 
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o. (3.11; 



Using the expression of ei oc and the equality p^a A d n+1 x = dx v A J{z) a a™, the previous 
conditions read 

§£ = «(•>?-«>?),., % = ^*-AW 

A* 

This proves the lemma. ■ 

Derivation in local coordinates. It is also instructive to carry out the variational 
principle directly in local coordinates. Using the local expressions of and e, see (2.13) 
and (3.7), we obtain 

6 e = (J(<£) - p%a A - rf))d n+1 x + p» A dy A A d n x„ 

so that, writing locally = (x^,ip A (x^),ip A (x^),a^(x^)), we get 

^6 e = (/«) - r A {< - lu))d n+1 x + rj^rdx" A d% 



ill TA\"fl I fl /J" ~ ' TA Qx V 

K)) d^x. 



di) A 
dx^ 

Applying the variational principle in local coordinates yields 

*/ U A (^+ii-*i)+iK))d n+1 x 



u 



8^ 



dx» 

A 



dx^ 



it - d n+l x + 6o$ (-dx* A J(z) a + ^d n+l x^ 

+ sr - *zur B - d n+l x + [ 4-sr A ^ A )d^ x . 



I 

Therefore, since the last term vanishes, we get the conditions 



Sat K ' dx» M ,At ' dx» 

which consistently recover the conditions (3.10) obtained by writing the condition (3.9) 
in local coordinates. 

The following Theorem summarizes the results obtained in this paragraph. 

Theorem 3.6 The following statements for a smooth section ip : U C X — > Il g are 

equivalent: 

(1) if; is a critical point of the covariant Clebsch variational principle (3.8). 

(2) ^*(iyf2 e ) = 0, for all nx,n g - v Grtical vector fields on II g . 

(3) tp satisfies the equations (3.10). 
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3.2.3 Hamiltonian formulation 

Recall that, under the hypothesis that £ is hyperregular, the critical curve q(t) of the 
classical (as opposed to covariant) Clebsch variational principle (2.2), is necessarily the 
projection to Q of the canonical Hamiltonian flow on T*Q associated to the optimal 
Hamiltonian H(a q ) = H(a q ,£*(a q )), see (2.7). This is a particular instance of a well- 
known result in optimal control theory based on the Pontryagin maximum principle, 
see e.g. [1]. As already mentioned earlier, for the Clebsch optimal control problem, the 
optimal Hamiltonian coincides with the collective Hamiltonian h o J associated to the 
Lie-Poisson Hamiltonian h obtained from the cost function Lagrangian £ by Legendre 
transformation. 

We shall show that the same situation arises in the covariant case, namely, that 
the stationarity condition obtained from the covariant Clebsch variational principle ad- 
mits, in the hyperregular canonical covariant Hamiltonian formulation with re- 
spect to a Hamiltonian section h : n — > J l Y* . This section is obtained from the 
section e by inserting the momentum map relation obtained from one of the station- 
arity condition. This Hamiltonian section naturally involves the Hamiltonian density 
h : L(L(TX, g), A n+1 X) — > A n+1 X associated to £ by Legendre transformation, i.e. 

HO = «d n+1 x - £«), < = j^ a ( 3 - 12 ) 

To make clear the distinction with the Hamiltonian density h, we call h the covariant 
Lie-Poisson Hamiltonian. This terminology is justified by the fact that, from (3.12), 
the covariant Euler-Poincare equations (2.21) for i are equivalent to the covariant Lie- 
Poisson equations for h given by 

div v + adlh v = 0. (3.13) 



It will be convenient to define the following map, which is yet another formulation 
of the covariant momentum map. We define 

j : n -> L(L(TX, g), A n+1 X), := dx» A J(z) a , (3.14) 

where z G J 1 !^* is such that fi(z) = u. Note that j verifies 

Kx,y a ) =j(u,a), for all a £ L(TX, g), 

where j is the bilinear bundle map defined in (3.5). 

Theorem 3.7 Suppose that the Lagrangian density £ is hyperregular, and define the 
associated Hamiltonian density h : L(L(TX, g), A n+1 ) — > A n+1 X by Legendre transfor- 
mation. 

Then, a section ijj — (uj , a) : U C X — > U g is a critical point of the Clebsch variational 
principle if and only if the section u is a solution of covariant Hamilton's equations 

io*(i v n h ) = 0, for all V £ X v (n), 
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i.e., in local coordinates 

w )» "337 = "ITTtM' 



y .4 

relative to the Hamiltonian section h : II — > J X Y* defined by 

h(w) = s 7 (w) - 7r^ y /i(j(a;)). (3.15) 

Proof. We first check that inserting one of the stationary condition, namely 4^ = 

dx^ A J{z) a = }(oj)%, in the section e defined in (3.6) yields the desired expression h in 
(3.15). This is the covariant analogue of the passage from the Pontryagin Hamiltonian 
H to the optimal Hamiltonian H, recalled in (2.7)-(2.8). We have 

e(w, a) = (s 7 (w) - j(u, a) + ir xx £(a), a) 

= (si(u)-ir* XtY (0(u;),a)-£(o-)),a) 

= S 7 (w) -7T XtY h(j(u)), 

which recovers the Hamiltonian section (3.15). Locally, it reads 

h(x»,y A , PA ) = ((x^y A , P %-h loc (x^y A ,p^) 

with h loc (x fl ,y A ,p'X) = h(j(x^, y A ,p A ) v a ) -Pa1£ = KPaW^ ~ pfotf, where d a denote a 
basis of the Lie algebra g. We now compute the covariant Hamilton equations. Denoting 
by v£ the variable of h, and using the notation := -^p- for the functional derivative 
of h, we have 



9hi oc dh A A A A A A 

dh toc dh 



y Md a ) A B - P^(7 m A ),b = p1K(d a ) A - ( 7 ;f), B ) = p^a A - tf) 

From this we obtain that the canonical covariant Hamilton equations 

du A dh ioc du'X dh loc 



B- 



dx^ du A y ~" dx^ du 

coincide with the first two equations in (3.10). To check complete equivalence with the 
stationarity condition (3.10), it remains to check that the definition we took for cr^ in 
the present proof coincides with the third equation for cr^ in (3.10). This follows by 
using the Legendre transformation relation between h and I: 

BP 

h(K) = KaZd n+1 x-£(a«), with 



H 1 



The Hamiltonian section h in (3.15) will be referred to as the covariant collective 
Hamiltonian, since the relation (3.15) is the extension of the relation H = ho J to the field 
theoretic context. Note that, in general, an Ehresmann connection is needed to obtain 
a globally defined Hamiltonian section h from the covariant Lie-Poisson Hamiltonian h. 
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3.2.4 Clebsch variables and the covariant Euler-Poincare equations 



Recall from §2.1 that the classical Clebsch variational principle reproduces the Euler- 
Poincare equations for the cost function Lagrangian £(£). We now prove a similar state- 
ment in the covariant case. To obtain this result, we need to assume that the open 
subset U on which the sections ip are defined is sufficiently small in order to allow the 
choice of the local trivial connection on U, which is simply the natural projection whose 
action in coordinates is (0, Va), i.e., the components 7^ = 0. 

Theorem 3.8 Let ip = (lu, a) : U C X — > Il be a critical point of the covariant Clebsch 
variational principle 

S [ ^*e e = o. 

Ju 

Then the section o : U C X — > L(TX, g) is a solution of the covariant Euler-Poincare 
equations 

div ^- + ad; ^ = 0. (3.16) 
da oa 

Proof. The proof follows from a computation in local coordinates. Let £ : U Cl->£| 
be an arbitrary smooth function. Using the stationarity conditions (3.10) and the fact 
that we chose the locally trivial connection 7^ = 0, we have 



X 



^AK^]j L d n+L x = -dx» A J(z) a [a^i\ 
=-(ad fT — ,£ 



which proves the result. In the fifth equality we used the local expression of the Jacobi- 
Lie bracket of vertical vector fields on 7. In the sixth equality, we used the formula 
[£,y,Vy}jl = ~[^v]y for a left Lie algebra representation (see e.g. §9.3 in [36]) together 
with the definition of the mapping J . ■ 



Canonical VS noncanonical structures. Recall that, in classical mechanics, Hamil- 
ton's equations are sometimes called noncanonical when their underlying Poisson struc- 
ture is not associated to a symplectic form. A lot of effort has been devoted to recast 
noncanonical Hamilton's equations into a canonical (i.e. symplectic) form by introducing 
auxiliary variables, like the well-known Clebsch variables in fluid mechanics, [15], [16], 
[37], [33]. They represent a first attempt to understand these equations as Hamiltonian 
system, before the introduction of Poisson manifolds. 

In the covariant case, and in order to emphasize the analogy with classical mechanics, 
we will use the same terminology to characterize the covariant Hamilton's equations 
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arising in this paper: the covariant Lie-Poisson equations (3.13) will be referred to as 
noncanonical as opposed to the canonical covariant Hamilton equations (2.17) associated 
to the canonical multisymplectic form. 



Curvature condition. It is important to mention that the covariant Clebsch vari- 
ational principle only yields the covariant Euler-Poincare equations and that the zero 
curvature condition (2.22) may not be satisfied in general. In the Lemma below we make 
this statement more precise. 

Lemma 3.9 Assume that Y is a trivial bundle, Y = X x Q and let if) = (ou, a) be a 
critical point of the covariant Clebsch variational principle. Then the following relation 
is satisfied 

(da-[a,a]) Q (q)=0 (3.17) 

as a skew-symmetric linear map T X X x T X X — > T q Q. 

This condition is weaker than the zero curvature condition (2.22). In the particular 
situation when the G-action on Q is free, then (3.17) and (2.22) are equivalent. 

When dimX = 1, then both conditions are always satisfied. 

Proof. Let V be a torsion free connection on Q and let (t, s) G M 2 — > x(s, t) G X be a 
smooth function. We have the relation 



TJs~ 



s=0 



d 
Iff 



q(x(s,t)) 



D' 



t=o 



Dt 



t=o 



d 



q(x(s,t)) = 0, 



(3.18) 



where D^/Ds denotes the covariant derivative associated to V. Using the notations 
x(s) = x(s,0), x(t) = x(0,t), Sx(t) := £\ s=Q x(s,t), Ax(s) = f t \ t=Q x(s,t), and the 
relation dq(x) = o~Q(q(x)), we compute 



D 



Ds 



s=0 



d_ 

dt 



q(x(s,t)) 



t=o 



Pi 

Ds 



dq(x(s)) ■ 5x(s) 



D' 



s=0 



Ds 



(a(x(s)) ■ Sx(s)) (q(x(s)) 



dq(x)-Ax 



(a(x)(Sx)) Q (q(x))+(- 



a(x(s)) • Sx(s)) (q(x)) 



s=0 



so that the relation (3.18) reads 







V^a^i , ^(^X^nO?^)) - V/ w n ff (a(x)(Ax))Jq(x)) 



d 

ds 



s=0 



a(x(s)) ■ 6x(s) - 1 



a(x(t)) ■ Ax(t) (q(x)) 



(a(x) ■ Ax) Q , (a(x) ■ 5x) Q ^ (q(x)) + (da(x) ■ (Ax, 5x)) Q (q(x)) 



JL 



[a, a] + da) Q (q(x)) ■ (Ax, Sx), 



where we used the property [£q, tjq] = —[£, t]]q for the infinitesimal generators of left ac- 
tion associated to Lie algebra elements £, 77 G Q, and a property of the exterior differential 
of a one-form. Since 5x, Ax G T X X are arbitrary, the result follows. ■ 
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Covariant Clebsch variables. From Theorem 3.7 we know that if (u, a) is a critical 
point of the covariant Clebsch variational principle, then oj is necessarily a solution of the 
canonical covariant Hamilton equations for the collective Hamiltonian section. Further- 
more, from Theorem 3.8 we know that o verifies the covariant Euler-Poincare equations 
and therefore, by Legendre transformation, v = j(u>) is solution of the noncanonical 
covariant Lie-Poisson equations (3.13). 

This means that u is a Clebsch variable for the covariant Lie-Poisson equation system, 
the relation between the canonical and noncanonical equations being given by the map 

u e n h-> e l(l(tx, q ), A n+1 x). 

Recall that this map is nothing else than a reformulation of the covariant momentum 
map associated to an action of the Lie group Q on the total space Y of the fiber bundle. 

Note finally that given a Lie group Q and a manifold X, any fiber bundle ttx,y '■ Y — > 
X on which Q acts (with trivial action on X) provide Clebsch variables for the covariant 
Lie-Poisson equations on L(TX,g), by means of the covariant momentum map. 

The situation is summarized in the next theorem. 

Theorem 3.10 (Clebsch variables for the covariant Lie-Poisson equations) 

Consider a manifold X and a Lie group Q. Let I : L(TX, g) — > A n+1 X be an hyperregular 
Lagrangian density and denote by h the associated Hamiltonian density. 

Let hx,y '■ Y —t X be a fiber bundle over X and suppose that Q acts on Y and induces 
the trivial action on X. Let h be the collective Hamiltonian section (3.15) defined locally 
with the help of the trivial connection, that is, we have 

h(x»,y A ,jf A ) = {x»,y A ,p%-h{){x»,y A , P » A ))). 

If u : U C X — >■ II is a solution of the (canonical) covariant Hamilton's equations 
for h, then j(w) G L(L(TX, g), A n+1 X) is solution of the covariant (noncanonical) Lie- 
Poisson equations for h. 

4 Examples 

In this section we apply the covariant Clebsch variational principle (VP) to various 
particular cases of special interest. First we consider the case given by a trivial fiber 
bundle Y = X x Q and a Lie group action of G on Q. Then we present the special cases 
when Q is a vector space on which G acts by representation and the case when Q is a 
group on which G acts by translation. This approach allows us to develop the covariant 
analog of several Clebsch variational formulations of the Euler-Poincare equations in 
mechanics ([21]), such as the symmetric formulation of rigid bodies, the coupled double 
bracket equations, the Clebsch variables for perfect fluids, and the singular solutions of 
EPDiff equations. The approach developed here applies to any equations admitting a 
covariant Euler-Poincare formulation, such as the G-strand equations for various choice 
of Lie groups G and Lagrangian £, [32]. We consider in details the cases G = Diff(M) 
and G = Diff„ ;(M), the associated covariant momentum maps and covariant Euler- 
Poincare equations (called covariant EPDiff and covariant EPDiff vo {). In the first case, 
our approach allows for the construction of singular solutions (filaments, sheets, etc.) of 
the equations. 
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4.1 Covariant Clebsch VP on trivial fiber bundles 

We now consider in detail the case when the bundle 7Tx,y '■ Y — > X is trivial. We 
will denote by Q the fiber so that we have Y = X x Q. Note that a smooth section 
: U C X — > Y is of the form (j){x) = (x, q(x)), where q : U C X — > Q. 

Recall that the fiber of the first jet bundle is an affine space. In the case of a 
trivial bundle it reads J l (X x Q)( x ,q) — id,T x x x L(T x X,T q Q) so it can be canonically 
identified with the vector space L(T x X,T q Q). This corresponds to the isomorphism 
j l (t>{x) i — y dq(x) for a section = (x,q(x)). 

The dual jet bundle and restricted dual jet bundle are given by 

J\X x Q)* = L(L(TX, TQ), A n+1 X) x A n+1 X and II = L(L(TX, TQ), A n+1 X) 
= (TX ® T*Q O A n+1 X) x A n+1 X = {TX ® T*Q ® A n+1 X) 

so that the line bundle /i : J X (X x Q)* — > II is trivial, and therefore Hamiltonian sections 
can be identified with maps h : II — > A n+l X covering the identity on X. 

Covariant momentum maps. Recall that in the field theoretic context, covariant 
momentum maps are defined on the dual jet bundle and take values in the bundle 
L(g,A n J\XxQY). 

In the special case of lifted actions we have seen in (2.18) that the value of the 
covariant momentum map is the pull-back of an expression in L(g,A n (X x Q)). This 
was the motivation for the definition of the map J : J 1 (X x QY — > L(g,A n (X x Q)) 
in (3.4). In the more special case of a group action on Y covering the identity on X, 
further simplifications arise in the expression of J . In view of the its occurrence in the 
stationarity condition of the Clebsch variational principle, it was natural to reformulate 
the covariant momentum map in this case as a map 

j : L(L(TX,TQ),A n+1 X) ^ L(L(TX,g),A n+1 X), 

see (3.14). In the case of a trivial bundle, the ^-action on Y = X x Q is the one naturally 
induced by a given ^-action $ on Q. We know check that the map j has a very simple 
expression in terms of the momentum map J : T*Q — > q* of the cotangent lifted action 
$ T * : G x T*Q T*Q of $. Indeed, using (3.4) and (3.14), we have 

(j(tu), a) = = dx» A J(z)M = dx» A P v A a^d n x v = P^d n+l x. 

So, for uj = v x <E> ctq ® fix and a x € L(TX, $j), we get 

(j(v x ® a q ® Hx), cr x ) = (a, a x (v x ) Q ) [i = (J(a 9 ), a x (v x )) \i x = (v x g> J(a 9 ) ® fx x , a x ) . 

Therefore, we have the formula 

j(v x <g> a q ® fix) =v x ® J(a 9 ) ® [i x . 

It will be also useful to write ) in coordinates as follows 

= ® c// ® rf n+1 x) = ® 3(dy A ) ® rf n+1 x 

= 9 M <8> J(w^/) ® rf n+1 x = ® J(w M ) (8) rf n+1 x. (4.1) 
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The covariant Clebsch variational principle. In the trivial bundle case, the vector 



bundle n g reads 



n g = L(L(TX, TQ), A n+1 X) x XxQ L(TX,g) 



and we will denote by ip = (u, a) : U C X — > H g the smooth sections of this bundle. 
Note that such section naturally encodes a smooth map q : U C X — > Q, given by the 
projection of u on Q. In other words, we have u(x) G L(L(T X X, T q r x }Q), A X +1 X). 

Using these notations, one easily verifies that the covariant Clebsch variational prin- 
ciple (3.8) reads 

5 [ e(a) + (u,dq-a Q (q)) = 0, 
Ju 

over the space of smooth sections if) — (u, a) : U C X — > Il . 

A detailed explanation of the formula appearing under the integral is in order. The 
symbol dq denotes the differential of the map q : U C X — > Q induced by uj, therefore 
dq(x) G L(T x X,T q ( x )Q). We have a(x) G L(T X X, q) and the notation dQ^q) denotes 
the function U C X — > L(TX,TQ) sending x G U to the linear map v x G T X X i— > 
( cr (i ; x))Q(9) ^ T q ^Q. The duality pairing denotes the natural duality between the 
spaces L(T x X,T q Q) and L(L(T X X, T q Q), A" +1 X). In local coordinates, we thus have 



6 [ t{o$ + ^ A {d,q A - a*)d n+1 x = 0. 
Ju 



A direct computation or an application of Lemma 3.5 yields the stationarity conditions 

Thanks to the expression (4.1), the last condition can be rewritten as 

S£ 



3(io^)d n+1 x. 



We now derive these conditions by using a global geometric formulation. To achieve 
this goal we need to fix a covariant derivative V on Q. For simplicity, we assume that 
V has no torsion. This covariant derivative induces a differential operator, denoted V, 
acting on smooth functions / : X — > TQ and defined by 



v 



/(*)), 



£ = 



where c{e) G X is a curve such that ^| £ _ C ( £: ) = v x and D v / De is the covariant 
derivative associated to V of the curve e (->■ f(c(e)) in TQ. Note that V Vx f\x) G T q ^Q, 
where q : X — > Q is the map induced by /. Given a variation ip £ = (u £ , a £ ) of the section 
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(uj, a), we compute 
d_ 

de 



e=0 



/ £(a) + {uj,dq-a Q (q)} 
Ju 



+ 



5a 
Ta M 

M A 

— , da 
da 



51 •( w 



e=0 



(^)q(&)} 



e=0 



+ ( W 



dq £ 



a F 



Me 



£=0 



+ (5 v u, dq - a Q (q)) + (u, V5q - (Sa) Q (q) - V Sq a Q 



5 v u, dq - a Q (q)) - (div V (u;) + (u, Va Q ) , 5q 



(4.3) 



In the second equality we used the general formula 



d 



D 



De 



-a F ,e F 



a F 



De 



where e(e) £ E is a curve in the total space of a given vector bundle E — > X endowed 
with the covariant derivative V B , and a(e) £ E* is a curve in the total space of the 
dual vector bundle E* — > X endowed with the induced covariant derivative V B *. The 
curves e(e) and a(e) are assumed to cover the same curve x(e) £ X. In the present case, 
the vector bundle is E = L(TX,TQ) — > X x Q and the covariant derivative is the one 
induced from the covariant derivative V on Q and a covariant derivative V x on X. One 
checks that our expressions do not depend on V x . For simplicity, we thus also denote 
by V the covariant derivatives on L(TX,TQ) and its dual. 

u £ and we use the following formulas. 



In the third equality we defined S v tu := yj- 



De 



e=0 



Pi 

~De~ 



Pi 
~De~ 



dq £ (x) ■ v x 



E = 



Pi 

De 

Pi 
Dt 



£=0 



(1 

dt 



Qe(x(t)) 



t=0 



Pi 

Dt 



t=o 



d_ 

de 



6 = 



Sq(x(t)) = V Vx 8q{x) 



t=0 

D v 



e=0 



e=0 



- Va ? (o-(z) ■ v x ) Q (q(x)) 
} {q)){x)-v x , 



De 

= {5a{x) ■ v x ) Q (q(x)) 
= (Nq(?) + V J9 (<t) 

where the last equality shows how one applies the section (Sa)Q(q) + V^cx),^) to v x . 

In the fourth equality, we used Stokes theorem together with the definition of the 
covariant divergence div associated to V and the fact that variations vanish at the 
boundary. We refer to [19] for more details about the operators V and div . 

From (4.3) we thus obtain the intrinsic formulation of the stationarity conditions 
(4.2): 

,r v, , , _ . . 6£ 



div" (w) = - (u, V<Jq) , dq = a Q (q), — = j(w). 



(4.4) 
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Clebsch variables and collective Hamiltonian. We first compute the collective 
Hamiltonian of the covariant Clebsch principle. Since the bundle is trivial, we can 
choose the trivial connection, i.e. the one given by the horizontal subbundle H (X x Q) — 
TX x {0} C T(X x Q). Recall also that the line bundle ft : J\X x Q)* ^ U is trivial, 
and therefore Hamiltonian section h : IT — > J X Y* = U XxA n+1 X can be globally written 
as h(a;) = (u, — H(u)), where H : IT — > A n+1 X is a smooth map covering the identity on 
X. With this identification, we obtain that the relation (3.15) is equivalently written as 

H(w) = h(j(u>)) G A n+1 X 

We recognize clearly here the generalization of the collective Hamiltonian relation 

H(a) = h(3(a)) 

to the covariant case. 

From the expression of j in the case of a trivial bundle, we obtain that the covariant 
Clebsch variables for the noncanonical covariant Lie-Poisson equations for h are provided 
by the variables uj\ through the mapping 

u = oj^ (8) dq A ® d n+1 x ^)(uj)=d fl ® J (of) ® d n+1 x. 

The variables evolve according to the canonical covariant Hamilton equations for h. 
For completeness, we present here a global formulation of these equations. A convenient 
way is to use the covariant Hamilton's phase space principle 

5 [ u*h*Q = 0. 
Ju 

In our case, since h(u) = (u, — H(cj)), we have u*h*Q = oj ■ dq — H(w), so that the 
covariant Hamilton's phase space principle reads 



5 [ wdg-H(w) = 0. 
Ju 



Similar computations as the ones used in (4.3) yield the equations 

dq = -, dw U = - — , (4.5) 

where a covariant derivative V has been fixed on Q. Here dH/du is the fiber derivative 
of H and <9 v H/<9g is the partial derivative of H relative to the connection V. 
When H(u;) = h(](tu)), we get 



du \5uJ dq \ \5uJ 



Qi 



where v = In this case (4.5) consistently recovers (4.4). 

The properties of the covariant Clebsch variational principle in the trivial bundle 
case are summarized in the following proposition. 
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Proposition 4.1 (Covariant Clebsch variational principle on trivial bundles) 

Let G be a Lie group acting on a manifold Q and let X be a manifold of dimension n + 1. 
Let £ : L(TX, g) — > A n+1 X be a Lagrangian density. Then we have the following state- 
ments. 

(1) The covariant Clebsch variational principle for a section ip = (u,a) : U Cl-lII, 
is 



5 



f £(a)-(u,dq-a Q (q)) = 0. 
Ju 



(2) The section if; = (u, a) is a critical section of the covariant Clebsch variational 
principle if and only if it verifies 

v S£ 
div (w) = - (u, Voq) , dq = a Q (q), —= (4.6) 

where V is a covariant derivative on Q with vanishing torsion. 

(3) The section a : U C X — > L(TX, g) is a solution of the covariant Euler-Poincare 
equations 

div ¥- + ad; -^ = (4.7) 
da da 

and verifies the compatibility condition (da — [cr, o"])g(g) = 0. 

(4) Suppose that I is hyperregular and define the associated Hamiltonian h : L(TX, g)* — > 
A n+1 X . Let H : II — > K' i+1 X be the collective covariant Hamiltonian defined by 
H = ho). Then (u>,a) is a solution of (4.6) if and only if u is a solution of the 
covariant canonical Hamilton's equations 

dq = — , div u = — . 4.8) 

ooj oq 

4.2 Covariant Clebsch VP for linear actions 

We now suppose that the Lie group G acts by representation on a vector space V. We 
denote by v t— > gv this action. The infinitesimal generator is denoted by £,v( v ) — £, v an d 
the cotangent bundle momentum J : T*V — > g*, is denoted by J(v,p) = v op. 

Let X be a (n + l)-dimensional manifold. Since T*V = V x V*, a section of the 
restricted dual jet bundle II — > X can be written as a couple (v,oj), where v : X — > V 
and bJ : X — > TX®V* ® A n+1 X. Writing u = d ' M <E> uj^ ® dx n+l , the covariant momentum 
map is given by j(v,u) = <9 M <S> v ® dx n+1 6 TX <g) g* ® A n+1 X. 

The covariant Clebsch variational principle reads 



5 / £(a) + (u), dv - av) = 



for maps ff:[/cl4 T*X®g, t;:[/cl^^,andu:[/cl^ TX ® 1/* ® A n+1 X. 
It yields the stationarity conditions 



5a 



)(v,oj), dv — av = 0, div uj — a^ = 
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which imply the covariant Euler-Poincare equations for £, together with the compatibility 
condition 

(do- - [a,a])v = 0. 

Assuming that the Lagrangian £ is hyperregular, with associate Hamiltonian h, the 
collective Hamiltonian is H(v,u) = h(d^ L ®vou^® d n+1 x) and the covariant canonical 
Hamilton's equations are 



dv = —, divw 

000 ov 



or, in terms of h, 



5h 5h .. . . ., . 

dv — — v, awu = - — u; , with u — j{v,ou). 
ov ovv- 

The case of G-strands. When X = M? 3 (t, s), we write a = ^dt+'jds, cu = (<9 i ®m + 
d s (g) n)dt Ads, v = (d t ® 7 + d s <E> X)dt A ds, with £(t, s),^(t, s) G Q, m(t, s),n(t, s) G V* , 
7(t, s), X(t, s) G g*. The covariant Clebsch VP becomes 



8 / 7) + ( m ) ^ ~~ € v ) dtds + (n, d s v — *yv) dtds = 0, 
Ju 

and yields the stationarity conditions 

— = (i> o m)ai A as, — = no njat A as 

04 07 

<9 t v = d s v = 7f , d t m + d s n = ^m + 777,, 

so that £, 7 verify the G-strand equations 

: + ad £ — r + <9 S — + ad — = 0. 

o£ s 04 07 ' ^7 

The compatibility condition is (dt'j — d s ^ — [£, j])v = and the Hamilton equations read 

a 5h a Sh a , a Sh 6h 
o t v = — v, o s v = —v, o t m + o s n = —m + -77-n. 
07 oa 07 oA 

Adjoint action and covariant coupled double bracket equations. A particularly 
interesting example is the case of the adjoint action, that is, we take V = Q and gm = 
Ad g m for g G G and m G g. Note that the infinitesimal generator is £g( m ) = ^d^m = 
[£,m], and the cotangent bundle momentum map is J(m,p) = mop = — ad^p. In 
this case, the classical Clebsch variational principle yields the coupled double bracket 
equations 

m = [m, [p, m]], p = [p, [p, m]] (4.9) 

of [4], as shown in [24], for an appropriate choice of the Lagrangian. We now present 
the covariant analogue of this system. 

Given a (n + l)-dimensional manifold X, the covariant Clebsch variational principle 
reads 

5 / £(a) + (u, dm - [a, m\) = 
Ju 
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for maps a : U C X ->■ T*X <g) g, m : U C X -> g, and w : C/ C X -> TX <g> g* g> 
A n+1 X, and where [a, m] G L(TX,g) is defined by [<7, m](v x ) := [cr(v x ), m]. It yields the 
stationarity conditions 

5£ 

— =)(m,u), d^m - [a M m] = 0, + ad^ = 0, (4.10) 

where j(m, cu) = — <9 M g) ad^a;^ £g> together with the compatibility condition 

[da — [a, a],m] = 0. 

Remarkably, the last two equations in (4.10) can be rewritten as 

d CT m = 0, div CT u; = 0, 

where we interpret a as a connection on the (left) trivial principal bundle IxG->G, 
and where d°" and div a denote, respectively, the covariant differential and covariant 
divergence associated to a . 

In the hyperregular case, these equations can be written with the help of the Hamil- 
tonian h as 

d 6h/5v m = 0, di\ Sh/Sv u = 0, where u = j(m,u). 

Suppose that g is endowed with an Ad-invariant inner product 7 and fix a Riemannian 
metric g on X. In this case, one can naturally identify the dual L(L(TX, g), A n+1 X) 
of L(TX,g) with L(TX,g) by using the duality pairing ((cr, 77)) := g^ v ^ a bV a jt v - The 
stationarity conditions are 

5£ 

— = [m, 00], dm — [a, m] = 0, div 97 00 — Tr 9 [<T, 00] = 0, 

where the functional derivative 81/ 5a G L(TX,g) is associated to the pairing ((,)), 
div 97 : T(L(TX, g)) — > ^(X, g) is the divergence associated to g and 7, and Tr 5 is the 
trace relative to g. If, moreover, the reduced Lagrangian is £(a) = ^\\a\\ 2 , where || || is 
the vector bundle norm associated to g and 7 on L(TX,g), then we have 

[m,w M ] = a 

so that, inserting this in the equations above, we get 

dm — [[m, u)],m] = 0, div 97 uj — Tr s [[m, u], uj\. 

These are the covariant coupled double bracket equations which extend (4.9) to the 
covariant case. They can also be written using covariant differential operators as 

d M m = 0, div^ 1 u = 0. 

Note that the associated covariant Euler-Poincare equations read 



div 57 ~ Tr s 



°'Ta 
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For the choice £(a) 
coupled double bracket equations imply 



\ || a || 2 , the second term vanishes and therefore the covariant 



div 97 a 



0. 



The results obtained in this paragraph, which extend those of [4] and [24] (see esp. 
equation (31)) to the covariant case, are summarized in the next proposition. 



Proposition 4.2 (Coupled double bracket formulation of covariant EP equations) 

Let i : L(TX,q) — > A n+1 X be a Lagrangian density assumed to be hyperregular, and 
consider the associated Hamiltonian density h. Assume that m : U C X — >• g, and 
oj : U C X — >• TX <S> g* <E> A n+1 X verify the canonical covariant Hamilton equations 



dm — 



5h 



S[m, ui] 



m 



0, diVgryU — Tl g 



5h 



8[m, oj] 



,oj 



Then a, given by = j(m, u), verifies the covariant Euler-Poincare equations, together 
with the compatibility condition [da — [a, o~], m] — 0. 



4.3 Covariant Clebsch VP for Lie group translations 

Recall that in the special case when a Lie group G acts on itself by left translation, 
the cotangent bundle momentum map is given by right translation at the identity: 
J : T*G ->■ fl*, 3(a g ) = a g g~\ [36]. 

The covariant analogue is obtained by taking Y = X x G and letting G act on itself 
by left translation. We have £ G {g) = TR g (£) = £g and j : U -> L(L(TX, g), A n+1 ) reads 

j(v x <g> a g ® ji x ) =v x ® a g g~ x ® /i x . 

A section oj of the restricted dual jet bundle LT — > X induces a map g : X G such 
that u>(x) E L(L(T x X,T g{x) G),A n+1 X). 

The covariant Clebsch variational principle reads 

5 [ i{a) + (u,dg-ag)=0, (4.11) 
Ju 

for sections u : U C X -»■ L(L(TX, TG), A n+1 X) and cr : U C X -> L(TX,g), and 
where dp(a;) G L{T x X,T g ^G) denotes the tangent map to g : X — > G and (crg)(x) G 
L(T x X,T g (x)G) is defined by (tr^)(i; x ) = (a(v x ))g(x) = TR g{x) (a(v x )). 

One can obtain the stationary conditions together with the Hamiltonian formulation 
by applying Proposition 4.1 and fixing a covariant derivative on G. We will however 
take the advantage of having a Lie group to trivialize the variational principle, by using 
the trivialized section v := ujg~ l : U C X — > L(L(TX, g), A n+1 X). Using this section, 
the variational principle (4.11) can be rewritten as 

5 [ £(a) + (u, (dg)g- 1 - a) = 0. 
Ju 
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A direct computation yields 



5 [ £(a) + (v, (dg)g- 1 - a) 
Ju 

and we get the stationarity conditions 
81 

— = v , dg-ag = 0, div(z/) + ad^-i 1/ = 0, 

where v = oog^ 1 . Note that from this expression of the stationarity conditions, it is clear 
that they imply the covariant Euler-Poincare equations, by simply inserting in the third 
equation, the expression for a and v given from the first two equations. 

Since the action is transitive, the compatibility condition coincides with the zero 
curvature condition (2.22), see Lemma 3.9. 

In the hyperregular case, the collective Hamiltonian is H(u;) = hiug" 1 ) and the 
process of passing from the covariant canonical Hamilton's equation for H to the covariant 
Lie-Poisson equations for h, coincides with the covariant Lie-Poisson reduction process, 
see [12]. 

The next proposition summarizes the properties obtained so far. 

Proposition 4.3 Let £ : L(TX,q) — > A n+1 X be a Lagrangian density and suppose that 
G acts on Q = G by left translation. Then the associated covariant Clebsch variational 
principle can be written as 

5 [ £(a) + {v, (dg)g- 1 -a)=0, (4.12) 
Ju 

where a: U C X -»■ L(TX,g), g : U d X ^ G , and v : U d X ^ L(L(TX, g), A n+1 X). 
The stationarity conditions are 

5£ 

— = u, dg-ag = 0, div(z/) + ad^-i u = 0, 

and imply the covariant Euler-Poincare equations for £ together with the zero curvature 
condition (2.22). 

Example: convective strand dynamics. When the base manifold is X = IR 2 3 
(t, s) and the Lie group is the special Euclidean group G = SE(3) 3 (A, r), we can write 
a = (u, i)dt + (fi, Y)ds and v = (d t <g> (m, n) + d s <g> (M, N))dt A ds, where 

(u, 7), (Q, T) : IR 2 — > se(3) and (m, n), (M, N) : M 2 ->■ se(3)*. 
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In this case the covariant Clebsch variational principle (4.12) reads 

5 J £(u, 7, ft, T) + <(m, n), (A, ^(ftA, d t r) - (w, 7 )> dtds 

+ <(M, iV), (A, r)- 1 ^, 8 s r) - (ft, T)) dtds = 0. 

For an appropriate choice of £, this principle yields the equations of motion for strand 
dynamics in convective representation ([38]), and recovers the variational principle de- 
veloped in [19]. The variable a = (u, r j)dt + (ft, F)ds involves the convective angular 
and linear velocities (u, 7) and the convective angular and linear strains (ft, T), whereas 
v — {pi (g) (m, n) + d s Cg> (M, N))dt A ds involves the corresponding convective momenta 
(m, n) and (M, A r ). We refer to [19] for the detailed derivation and alternative variational 
formulations of convective strand dynamics. 



Example: the covariant version of the symmetric rigid body equations. The 

so called symmetric rigid body equations 

Q = QU, P = PU, Ueso(N), (4.13) 

are an alternative formulation of N- dimensional rigid body dynamics, associated to an 
optimal control problem [4], [6], [24]. Given an hyperregular Lagrangian i : so(A r ) — > R, 
for example the N- dimensional rigid body Lagrangian, one observes that if U is defined 
by 5£/5U = \{Q T P — P T Q), then U is a solution of the Euler-Poincare equations for £ 
on so(N). We shall now present the covariant analogue of this statement. 

Consider the right action of SO(N) on the trivial bundle X x GL(N) — > X by 
translation on GL(N). The covariant Clebsch variational principle is 



5 [ £{v) + (to, dQ - Qv) = 0, 
Ju 



for sections u : U C X L(L(TX,TGL(N)), A n+1 X) and v : U C X ->■ L(TX,so(N)), 
and yields the stationarity conditions 

= j(Q,^ = ^(Q^-KfQ), dQ = Qt;, diva; = a;t;. 
ov/j 2 

The last two equations are the covariant analogues of the symmetric equations (4.13) 
and imply the covariant Euler-Poincare equations on so(N) 



, . 5£ 
div — + 

ov 



0. 



This applies in particular to the 5'0(A r )-strand equations if X = IR 2 3 (t, s). In this 
case, using the notations v = Udt + Vds, and uj = (d t <8> M + 9 S (g) Ar)dt A ds, with 
C/(t, s), V(t, s) G so(iV) and M(t, s), N(t, s) e T Q ^ S )GL(N), the stationarity conditions 
are 

^ = J(Q T M - M T Q)dt Ads, 4? = 1(Q Tn ~ N T Q)dt A ds 
oU 2 ov 2 
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d t Q = QU, d s Q = QV, d t M + d s N = MU + NV 
and imply the SO (N)-stiand equations 



Si 



Si 



dt— + d s — + 



su 



sv 



'SU 



+ 



0. 



The particular case N = 3 is relevant for chiral models, in which case the Lagrangian 
density is i(U, V) = [\\U\ 2 — ~\V\ 2 )dt A ds, see [32] and references therein. 



4.4 The covariant EPDifF equations and singular solutions 

The covariant EPDifF equations. These equations are, by definition, the covariant 
Euler-Poincare equations in the special case when the Lie group is the diffeomorphism 
group Diff(M) of some manifold M. The terminology that we use comes from the same 
denomination (EPDiff) for the Euler-Poincare equations on the diffeomorphism group, 
[34]. The introduction of this class of partial differential equations, the covariant EPDiff 
equations, was motivated by the recent paper [32]. 

Given a (n + l)-dimensional manifold X we get from the general expression (2.21) 
of the covariant Euler-Poincare equations, the covariant EPDiff equations 

Si Si 

div x — + ad; — = 0, (4.14) 
ocr M oer^ 

where ad* m = Vm • u + Vu T ■ m + mdivu is the infinitesimal coadjoint operator 
associated to the Lie algebra g = X(M) of vector fields on M, where we identified q* 
with q using the L 2 -duality pairing. In (4.14), divx denotes the divergence of the X(M)*- 
valued vector field density £ G TX <g> X(M)* g> A n+1 X. Note that in these equations, 
a is a smooth section er : U C X — > L(TX, X(M)). This implies that, given iGl and 
m G M, we have <r(x,m) := cr(x)(m) G L(T x X,T m M). The zero curvature conditions 
read 

d x <r =[*,*], (4.15) 

where we used the notation to make clear that it means the exterior derivative relative 
to the variable x, that is, the exterior derivative of the one-form cr G fi 1 (X, X(M)). Note 
that (4.15) is an equality in fi 2 (X,X(M)). 

In the case X = IR 2 3 (t, s) and M = R, the covariant EPDiff equations recover 
the Diff (IR)-strand equations studied in [32]. Using the notation a(t,s) = u(t,s)dt + 
7(s, t)ds G j£(M), we get from (4.14) the Diff (R)-strand equations 

n Si , fl Si 51 51 

o t — + o s — + ad^ — + ad — = 0, 
ov 07 ov ' 07 

where ad* v = v'u + 2vu' is the infinitesimal coadjoint operator on £(R)*. The zero 
curvature condition (4.15) reduces to 

d t j -d s v= [u, 7] = ~fd m v - udmry. 



35 



In [34] a pair of momentum maps has been defined in the context of the (classical) 
EPDiff equations. These momentum maps are associated to the two natural diffeomor- 
phism actions on the space Emb(S', M) of embeddings of a k- dimensional manifold S into 
M. The left action of the diffeomorphism group Diff(M) yields the cotangent bundle 
momentum map Jl '■ T* Emb(5, M) — > X(M)* given by 

J L (Q,P)= / P(s)6(m-Q(s))d k s. 
Js 

This momentum map provides the expression for singular solutions of the EPDiff equa- 
tions. They may, for example, be supported on sets of points (vector peakons, dim S = 
0), one-dimensional filaments (strings, dim S = 1), or two-dimensional surfaces (sheets, 
dimS* = 2) when dimM = 3. The right action of the diffeomorphism group Diff (5) on 
Emb(5*, M) yields the cotangent momentum map 

3 R (Q,P) = P-VQeX(S)* 

which provides Clebsch variables for the EPDiff equations on X(S). Two Clebsch vari- 
ational principles can be naturally associated to these actions, [24]. 

We shall now present the covariant analogue of this setting and obtain singular 
solutions of the covariant EPDiff equations via covariant momentum maps. 



Covariant Clebsch variational principle I: singular solutions. We take the fiber 
bundle Y = X x Emb(5', M) and the principal bundle P = X x Diff (M), and consider 
the left action of Diff (M) on Emb(S', M) by composition. The restricted dual jet bundle 
is II = TX ® T* Emb(5', M) C*D A n+1 X. Using our general formula (4.1), we obtain the 
covariant momentum map )l '■ n — > TX ® X(M)* Cg> A n+1 X, 

j L (Q, w) = ® J(Q, u») ® d n+l x = ® [ ^(s)5{m - Q(s))d k s <g> d n+1 x. 

Js 

From Proposition 4.1, the covariant Clebsch variational principle reads 

5 [ + («, d x Q - o- Emb (Q)) = 0, (4.16) 
Ju 

over sections a : U C X -> L(TX,X(M)) and u : U C U ->■ II. As before, 
Q : U C X — > Emb(5', M) is the map induced by u}, that is, we have u>(x) G 
L(L(T x X,T Q(x) Emb(S,M)),A n+1 X). We will use the notations Q(x, s) := Q(x)(s) and 
oj(x, s) := uj(x)(s). In the variational principle (4.16), dxQ denotes the derivative of the 
map Q:UcX^ Emb(5', M), and the expression <r Emb (Q) G L(TX,TEmb(S, M)) is 
given by <r Emb (Q)(t; :r )(s) = (cr(v x ))(Q(x, s)) G Tq( XjS )M. The variational principle can 
be written as 

5 [ £(<t(x)) + (^(i),9,Q(i) - (^(x)) Emb (Q(x))} = 
Ju 

or, writing explicitly the duality pairing between TEmb(S,M) and its dual, 



'(x))+ / (uj"(x,s),d lJt Q(x,s) - o- M (x,Q(x,s)))c/ fc s 
Js 



0. (4.17) 
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The stationarity conditions read 
S£ 



[ u»{s)5{m - Q(s))d k s <9 M Q(x, s) = a^x, Q(x, s)) 
Js 

div x u: = -((V<t/oQ) • u/\ 



In the hyperregular case, eliminating the variable cr, these equations are equivalent 
to canonical covariant Hamilton's equations for the collective Hamiltonian H(Q, u;) = 
/i0i(Q, a;)) and read 

dQ = ^, drv^ = - — 

on the restricted dual jet bundle II. Of course, one has to assume that the Hamiltonian 
h converges on the value of the covariant momentum map ji, that is, the collective 
Hamiltonian H = h o ) L is well-defined. 

In this case, it follows from Proposition 4.1 that cr is a solution of (4.15) and satisfy 
the compatibility conditions (d^cr — [cr, cr]) Emb (Q) = 0, which can be written more 
explicitly 

(d x cr(u x , v x ) + [<r(u x ), cr(v x )}) o Q(x) = 0, for all u x , v x G T X X, (4.18) 

as an equality in Tqm Emb(5', M). 

The covariant Clebsch variational principle generalizes easily to the case when the 
bundle is 

Y = X x Emb(S,M) N 3 (x,Q a ), a = l,...,iV, 
in which case the covariant momentum map is 

N 

j L (Q a , u a ) = Y^d,® ^(s)8(m - Q a (s))d k s ® d n+1 x. (4.19) 
«=i Js 

In the Proposition below, we consider the case of a H 1 Lagrangian density, which is 
the natural generalization of the higher dimensional Camassa-Holm equations to the 
covariant case. 

Proposition 4.4 (Singular momentum map solutions of H 1 covariant EPDiff) 

Consider the Lagrangian density I : L(TX,X(M)) — > A n+1 X given by 



f^lkolllrx-iEW^)^ 



where ||u||^i = ||u|| 2 2 + a 2 ||Vu|| 2 2; a > 0. Then the corresponding covariant EPDiff 
equation is 

dpif + Vi/ 1 • cr^ + Vcr I ■ + v" div cr At = 0, 

with is — (1 — a 2 A)cr and v % = — (1 — a 2 A)crj ; i = 1, n. 
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This equation admits the singular solutions cr = cr^dx^ given by 

N 

cr (x,m) = V] / s)G(m, Q a (x, s))d k s 

«=i ^ s 

N 

(Ti(x,m) = -V / u/ a (x,s)G(m,Q a (a;,s))d fc s, z = l,...,n, 
a=i ^ 

where G is the Green's function for the linear operator (1 — a 2 A) and where (Q a , w a ), a = 
1, ...,N verify the canonical covariant Hamilton's equations 

dQ"=g^, div x a, . = -_ , „=!,...,*, (4.20) 

/or i/ie collective Hamiltonian 

N N 

H((Q a , a; a )f =1 ) = ^ E 2 // s)G (Qa(x ' s) ' Q&(x ' S ' )} S ' )d " S rfV 

^=0 a,b=l J •' S 

on the restricted dual jet bundle U = TX ®T* Emb(S r , M)^ <g> A n+1 X. 

Proof. These results follows from the discussion before the proposition. In particular, 
the expression for <x and cr^ follows by inverting the operator (1 — a 2 A) and using the 
associated Green's function. The expression for the collective Hamiltonian is obtained 
by using the formula H = h o ) L and noting that the covariant Lie-Poisson Hamiltonian 
h associated to I reads 

= - / J v fl (m)G(m,n)v^(n)dm dn. ■ 
h=o 2JJm 



Remark 4.5 (1) Note that any covariant EPDiff equations whose Lagrangian density £ 
is such that the collective Hamiltonian H = hoj L is well-defined, admits these momentum 
map singular solutions. 

(2) Note also that these singular solutions do not satisfy the zero curvature condition 
(4.15) but only the compatibility conditions (4.18), see Proposition 4.1. Therefore, 
one cannot reconstruct the dynamics of the sections x 6 U C X i— > (x, <fi(x)) G X X 
Diff(M) of the bundle X x Diff(M) — > X. The dynamics of these singular solutions 
takes naturally place on the restricted dual jet bundle n = TX <g> T* Emb(S f , M) u ® 
A n+1 X 3 (Q a ,u> a ) and is given by the canonical covariant Hamilton equations for H. 
This remark holds in general for covariant Euler-Poincare equations obtained through 
covariant Clebsch variational principles, see Proposition 4.1. 

(3) When dim M — 3, for example, the components of the singular solutions may be 
supported on sets of points (vector peakons, dimS 1 = 0), one-dimensional filaments 
(strings, dimS* = 1), or two-dimensional surfaces (sheets, dimS* = 2). 
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The case of Diff (M)-strands. When X = M 2 3 (t,s) and M = R, we recover the 
results of [32]. In this case, the manifold Emb(S f , M) is replaced by M. (the manifold S 
is reduced to a point). We use the notations 

<r = vdt + "yds and u a = (d t (§> M a + d s ® N a )dt A ds, 

where <r(i, s), 7(i, s) G X(R) and M a (t, s), N a (t, s) G Tg^^R. The covariant momen- 
tum map (4.19) is therefore given by 

TV 

j L ((Q a , u;X=i) = (9t ® M„ «J(m - Q a ) + d s <g> iV a <f(m - Q a )) dt A ds 

a=l 

and recovers the peakon solutions 

N N 

u(t, s,m) = Yl M - G ( m > Q a (t, s)), 7 (t, s, m) = - ^ iV a G(m, Q a (t, s)) 

a=l a=l 

of the Diff (M)-strand equations found in [32]. In particular, our approach yields a natural 
interpretation of the equations governing the dynamics of these singular solutions. In- 
deed, from (4.20), the fields Q a , N a , M a obey the canonical covariant Hamilton equations 
for the collective Hamiltonian H given by 

1 N 1 N 

H((Q a ,iV a ,MX =1 ) = - N a G(Q a , Q b )N b + - ^ M a G(Q a ,Q b )M b 

a, 6=1 a, 6=1 

on the restricted dual jet bundle II = TR 2 (g) T*R ® A 2 R 2 . These equations are 
d t Q a {t, s) = u(t, s, Q a (t, s)), <9 s Q a (£, s) = 7 (t, s, Q a (t, s)), 

«9G 



<9 4 M Q + «9 s iV a = - ^(iV^ + M a M b ) — (Q a , Q b ). 

6=1 ^ 

They coincide with the equations derived in [32], interpreted there as a constrained 
classical Hamiltonian system, rather than an unconstrained covariant Hamiltonian sys- 
tem. 

The compatibility conditions (4.18) reduce to (d t j — d s v — "fd m v + vd m ~{) (Q a (t, s)) - 
0, for all a = 1, N and yield 

N 

J2(dtN b + d s M b )G(Q a ,Q b ) 
6=1 

N ( dG dG \ 

+ J^(M b N c -N b M c ) i^(Q a ,Qb)G(Q a ,Q c ) + -^- b (Q a ,Q b )G(Q b ,Q c )j =0, 

for a = 1, N. The same condition can be derived by writing d ts Q a = d st Q a . 
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For completeness, we rewrite now the covariant Clebsch variational principle (4.17) 
in the particular case of Diff (M)-strands. It reads 



5 [ £(u(t, s), 7 (t, s)) + Y] M a (t, s){d t Q a {t, s) - u(t, s, Q a {t, s))dtds 

N 

+ ^ N a {t, s)(d s Q a (t, s) - 7 (t, s, Q a (t, s))dtds = 0. 



a=i 



Covariant Clebsch variational principle II: covariant Clebsch variables. In 

order to implement the second covariant variational principle, we take the same fiber 
bundle Y = X x Emb(5, M) as before, and consider the principal bundle P = X x 
Diff (S), where Diff (S) acts on Emb(S', M) by composition on the right. From our general 
formula (4.1), we obtain the covariant momentum map )r : II — > TX (g) X(S)* <8> A n+1 X, 

j fl (Q, u) = d l ,® J(Q, ^) <g> d n+1 x = <8> u}» ■ VQ ® 

The covariant Clebsch variational principle reads 

5 / £(<r) + (u>, d x Q + VQ • cr) = 0, (4.21) 

over sections cr : [/ c X L(TX,X(S)) and w : C/ C X ->■ n, where Q : C/ C 
X — > Emb(S', M) is the map induced by uj. In the variational principle, VQ denote 
the derivative of Q(x, s) relative to s, and VQ • cr e L(TX, T Emb(S f , M)) is given by 
(VQ • cr)(t> x .) = VQ ■ (cr(ur)) G Tq( x ) Emb(S', M). The variational principle can be 
written as 

5 [ i{a{x)) + (w"(aO, 9 M Q(x) + VQ • cr M (x))} = 



or, more explicitly, 



5 [ £((t(x)) + [ (u^x, s), 9 M Q(x, s) + VQ(x, s) ■ <t m (z, s))) rf fc s 



o. 



As above, in the hyperregular case the critical sections (Q, u>) of this principle satisfy the 
covariant Hamilton equations for the collective Hamiltonian H = hoj R and cr satisfies the 
covariant Euler-Poincare equations for i on L(TX,X(S)), as well as the compatibility 
conditions (3.17) given here by 

VQ-(d x cr-[cr,cr]) = 0. 

The detailed computations are left to the reader and are obtained by extending those 
of [24, §4] to the covariant case. The sections (Q, u>) are interpreted here as covariant 
Clebsch variables for the covariant EPDiff equations, through the relation 

^- = -j fl (Q, w) = -d„ ®w"-VQ® d n+1 x. 

OCT 

Note that (4.21) differs from the variational principle of Proposition 4.1 by a change 
of sign in front of the infinitesimal generator. This change of sign is due to the use of a 
right action whereas Proposition 4.1 is formulated for left actions. 
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A pair of covariant momentum maps. The left and right actions of the diffeomor- 
phism groups on X x Emb (S,M) induce a pair of covariant momentum maps on the 
restricted dual jet bundle IT: 



n 




TX ® 3L{M)* <g> A n+1 X TX ® X(S)* ® A n+1 X, 

similar to the dual pair of momentum map considered in [3 1] and further studied in [25], 
the left leg being associated to singular solutions and the right leg to covariant Clebsch 
variables. 



4.5 The covariant EPDifF W0 / equations 

The covariant EPDiff„ / equations. We now consider the incompressible version of 
the covariant EPDiff equation studied above, that is, the case of the Lie group Diff^(M) 
of volume preserving diffeomorphisms of a Riemannian manifold M (assumed to be 
compact and without boundary, for simplicity). The Lie algebra of this group is given 
by the space X vo i(M) of divergence free vector fields on M. Making the identification 
X vo i(M)* = X vo i(M) and using the L 2 pairing associated to the Riemannian metric 
on M, we have the expression ad* m = P(V u m + Vu T ■ m), for all u, m e X vo i(M), 
where P : X(M) — > X vo i(M) denotes the orthogonal Hodge projector onto divergence 
free vector fields. 

Given a (n + l)-dimensional manifold X we get from (2.21), the following covariant 
Euler-Poincare equations on Diff^(M): 

5£ 5£ j 5£ 

dwx j- + V CTm1 — + Vcr r ■ — = gradp, (4.22) 
da oa^ p ba^ 

where a is a smooth section a : U C X -»■ L(TX,X(M)) and p : X x M ->• E. The 
symbol V and grad denotes, respectively, the Levi-Civita covariant derivative and the 
gradient on M, both computed relative to the Riemannian metric on M. If the first 
cohomology H 1 (M) of M is trivial, we can alternatively identify the dual Lie algebra 
with the space of exact two-forms on M, that is, X vo i(M)* = VL 2 ex (M). In this case, the 
equations take the vorticity form 

51 5£ 

div x — + £ — = Q, (4.23) 
da 5a ^ 

where £ a ^5£/5a^ denotes the Lie derivative of Si/Sa^ eTX® VL 2 ex (M) ® A n+1 X with 
respect to M, along the vector field a fl . 

In the 2D case, on can define a covariant stream section ip : U C X — > L(TX, jF(M)) 
for a and rewrite the equations as 



div ^ + fe^ = ' (4 - 24) 
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where { , } denotes the Poisson bracket associated with the volume form on M, viewed 
symplectic form. 



Covariant Clebsch variational principle and Clebsch variables. Recall that 
the classical Clebsch variables for perfect fluids arise from the momentum map J : 
T*F(M) -»■ X(M)* = £l 2 ex (M), J(A,/i) = dA A d/i associated to the cotangent lifted 
action the right composition by Diff(Af) on J-"(M), [37]. 

In the covariant case, we take the bundle Y = X x .F(M), and we get the restricted 
dual jet bundle II = TX ® T*J r (M) eg) A n+1 X 3 (A, u). The covariant momentum map 
is j(A, u) = £g> dA A dw M £g> c? n+1 x and the covariant Clebsch variational principle reads 

'(a) + (to, d x X + dA ■ a) = 0, 

over sections a : U C X ^ L(TX,X vol (M)) and (A, w) : U C [/ -»• II, where d x A 
denotes the derivative of the map A : U C X — > J-~(M). The stationarity conditions are 

51 

— =)(\,u), d x A + dA-cr = 0, div x u + Tr(dw • a) = 0. 

OCT 

In the hyperregular case, the first relation can be written as tr = gJ£^\ ■ The last two 
equations can be written as a covariant canonical Hamiltonian system 

= q^i dw x w = 

for the collective Hamiltonian H = /i o j on II, and yields the covariant EPDiff„ / equa- 
tions (4.23), together with the compatibility condition dA • (d^cr — [c,(t]) = 0. The 
sections (A, u) are thus interpreted as covariant Clebsch variables since they obey canon- 
ical covariant Hamilton's equations and yields the noncanonical covariant Lie-Poisson 
equations for h, in a similar way with the classical Clebsch variables (A,/i). 

For example, with the L 2 Lagrangian density £(cr) = Y^=o ^II^mIIl 2 we nave = 
der^ and the covariant Hamilton's equations for Clebsch variables become 

d x A + dA ■ d~ x (dA Aw) =0, div x u + Tr{du ■ d~\d\ A u)) = 0, 

where d _1 denotes the inverse of the exterior differential seen as an isomorphism d : 
Ql ol (M) —> VL 2 X (M). In the two-dimensional case, see (4.24), these equations are 

d^X + {A, {A, = (no sum on /i), d^u^ + {u^jA, uj^}} = (sum on /i) 

and yield a covariant coupled double bracket formulation of the two-dimensional EPDiff„ ; 
equations. They are reminiscent of the coupled double bracket formulation obtained in 
Proposition 4.2 in the context of adjoint actions, and extend to the covariant case the 
coupled double bracket formulation of perfect 2D fluids obtained in [24, §6]. 
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5 The covariant Pontryagin variational principle 

In this section we develop a general covariant variational principle that generalizes the 
covariant Clebsch variational principle and reproduces, as general case, the classical 
and covariant Hamilton-Pontryagin variational principles of [39] and [40], respectively, 
together with variational principles arising in geometric optimal control problems, e.g. 
[29]. 

Given a fiber bundle ttx,y '■ Y — > X, we consider another fiber bundle tty,b : B Y 
over Y with local coordinates y A , b a ), and define the bundles 

3 : = J 1 ^ x Y B and % = Tlx Y B ^Y, (5.1) 

referred to as the generalized Pontryagin bundle and restricted generalized Pontryagin 
bundle, respectively. 

The line bundle \i : J X Y* — > U induces naturally a line bundle m : 3 ~ * locally 
given by (x M , y A , 7r , p^, b a ) — > (x^ 1 , y A ,p^, b a ). The canonical forms 6 G r2 n+1 (J 1 F*) and 
n G Q^i^Y*) induce the forms 6 3 := vr} ly , 3 G fi n+1 (3) and fi 3 := 7r} ly , 5 Q G 

The variational principle that we define is associated to a smooth section e : — > 3 
called a generalized energy section, written locally as 

(x»,y A ,p%b a ) i y (x»,y A ,-e loc (x»,y A ,p%b a ),p%b a ). 

By definition, a generalized energy density is a smooth map £ : 3 — A n+1 X covering 
the identity on X and such that 

ix(d£ + fi 3 ) = 0, for all X G 3^(3), (5.2) 

where £ y (3) = {X G £(3) | ^nx o X = 0} is the space of all m-vertical vector fields on 
the generalized covariant Pontryagin bundle. From the condition (5.2) we obtain that a 
generalized energy density is locally given by 

£(x»,y A ,p%n,b a ) = (n + e loc (x»,y A ,p%b a ))d n+1 x. 

This definition of generalized energy density is the natural extension of the definition 
(2.15) of Hamiltonian densities. Similarly to this case, there is a bijective correspondence 
between generalized energy densities and generalized energy sections. This correspon- 
dence being given by 

ime = £ -1 (0). 

The Cartan (n + l)-form and (n + 2) forms associated to a generalized energy section 
e are defined by 

9 e := e*6 3 G Q n+1 {¥) and Q e := e*Q 5 G Q n+2 {^) 
and we have the relations 

m*Q e = Q z — £ and m*f2 e = Vt z + d£. 
One verifies the first equality by checking that both terms are locally given by 

-e {oc (z", y A , p% b a )d n+1 x + p\dy A A d%. 
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Definition 5.1 (The covariant Pontryagin variational principle) Letnxy '■ Y — > 
X be a fiber bundle and fix a fiber bundle hy,b '■ B — >■ Y and a generalized energy sec- 
tion e : *P — > 3- Let U C X be an open subset whose closure U is compact, and let 
if) : U C X *p be a local smooth section of nx,<$ '■ ^3 — > X . 
The covariant Pontryagin variational principle is 

5 [ -0*0 e = 0, (5.3) 
Ju 

for all variations ip e : U C X — > of ip (among smooth sections) such that ipo = ip and 
4>e\du — 4>\du, where <p £ is the section of ttx,y '■ Y — > X induced from ip £ . 

Locally, the variational principle (5.3) reads 

S I (^TT ~ eioc(x»,i> A ,r A ,r)) d n+1 x = 0, 



where we denoted locally the section as ip{x^) = (x M , ip A (x fJl ), ip^x^), ip a (x^)) 

Remark 5.2 (Formulation on 3) The variational principle (5.3) can also be formu- 
lated on sections ip : U C X — > 3 of the generalized Pontryagin bundle as 



s [ $*{e z -e) = o, 

Ju 



'U 

since we have ip*(Bz — £) = 4>*m*Q e = (tnot/;)*9 e = ip*Q e , where we defined the section 
ip := m o ip of Of course, the associated stationarity conditions only involve the 
projection ip :— m o ip of ip and the section ip is reconstructed from ip by imposing the 
zero energy constraint 

8 o$ = 0. 

In coordinates, the section ip is thus given by 

In the same way as in Theorem 3.6, we can prove the following result. 

Theorem 5.3 The following statements for a smooth section ip : U C X — > are 
equivalent: 

(1) ip is a critical point of the covariant Pontryagin principle (5.3). 

(2) ip*(ivQ e ) = 0, for all Ttxffl -vertical vector fields on *}3. 

(3) ip satisfies the following equations in local coordinates 



dip A dei oc dip^ A dei oc de 



loc 



dx^ dip A ' dxi 1 dip A ' dip a 



0. (5.4) 



We now show that this covariant principle unifies several known variational principles 
in field theories and classical mechanics, including the covariant Clebsch variational 
principle developed in this paper and the Hamilton-Pontryagin principle. 
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Recovering the covariant Hamilton-Pontryagin principle. Let C : J l Y — > 

A n+l X be a Lagrangian density and consider the fiber bundle B := J l Y — > Y . The gen- 
eralized Pontryagin bundle is thus given by = J Y* Xy J l Y . Define the generalized 
energy density S : J X Y* x Y J X Y A n+1 X by 

£(z, 7 ) = (z,j) - £(7). 

Since it is locally given by S(x tJ ', y A ,p^, v A ) = (jr + Pa v u ~~ ^{x^, y A , v A )) d n+1 x, it ver- 
ifies the condition (5.2). In this case, (5.3) recovers the covariant Hamilton-Pontryagin 
variational principle developed in [40]. Since e/ oc = p^v A — L(x^ , y A , v A ), the stationarity 
conditions (5.4) read 

^ = ,a d^ = dL_ ^_dL_ = Q 
dx» ^' dx» d^ A " dt/tf 

and consistently recover the Euler-Lagrange equations in implicit form. 

Recovering the covariant Clebsch variational principle. Let Q C Aut(Y) be a 
Lie group acting on Y by bundle automorphisms covering the identity. Let J : J l Y* ->■ 
L(g, A n Y) be the covariant momentum map and consider the bundle B := L(TX,q) — > 
Y. Given an Ehresmann connection 7, we consider the energy section defined by 

e(w, &) := (s 7 (w) - j(u, a) + ir* X Y l(a), a) , 

see (3.6). Then the covariant Pontryagin principle (5.3) recovers the covariant Clebsch 
variational principle of Definition 3.3. Since e/ oc = v\{o A — J A ) — K a a)} one verifies that 
the stationarity conditions (5.4) agree with (3.10). 

Recovering the covariant Hamilton phase space principle. Given a Hamiltonian 
density H : J X Y* — > A n+1 X, we consider the case when the bundle B — > Y is absent, 
i.e., 3 = J l Y*, and we take £ := %, so that e = h. In this case the covariant Pontryagin 
principle (5.3) recovers the covariant Hamilton phase space principle recalled in (2.16). 
The stationarity conditions (5.4) yield the covariant Hamilton equations (2.17). 

Recovering classical Pontryagin principles in optimal control and mechanics. 

In order to particularize the variational principle (5.3) to classical mechanics, we choose 
X := E, Y := R x Q -> E, and the bundle B := R x C -> E x Q, (t, c q ) ^ (t,q), 
where C — > Q is a fiber bundle over Q. In this case, the generalized Pontryagin bundle 
and restricted Pontryagin bundle defined in (5.1) are 3 = T*(R x Q) X(k x q) x(lxC) 
and *P = (T*Q x R) X( Rx q) x(E x C). A generalized energy section reads e(a q ,t,c q ) = 
(a q ,t, —e(t,a q ,c q ),c q ) and a local section of *P — > X reads ip{t) = (a(t),t,c(t)). With 
these notations, the variational principle (5.3) becomes 

5 I ?p*e e = 5 [ ((a(t),q(t)) - e(t,a(t),c(t)))dt = 0. (5.5) 
Jo Jo 

By choosing e(t, a q , c q ) := (a g , T(c g )) — L(c q ), where L : C — > R and T is a fiber bundle 
map r : C — > TQ covering the identity on Q, (5.5) recovers the variational principle 
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associated to geometric nonlinear optimal control problems via the Pontryagin maximum 
principle, see e.g. [29]. In this case C is the control bundle, L is the cost function of the 
problem, e is the Pontryagin Hamiltonian, and q(t) = T(c(t)) is the differential equation 
constraint. 

When C = TQ and e(t, a q , v q ) = (a q , v q ) — L(v q ), (i.e., T = icItq) then (5.5) recovers 
the Hamilton-Pontryagin variational principle ([39]); when C = Q x g and e(t,a q ,C.) = 
(o^'Cq^)) — ^(£)> 0- e - T(g,£) = £q (<?)), then (5.5) recovers the Clebsch variational 
principle ([24]); finally, when C is absent and e(t,a q ) = H(a q ), then (5.5) recovers the 
classical Hamilton phase space variational principle. 
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